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ABSTRACT 

A  pulse-echo  system  using  an  electrostatic  transducer  and 
steel  reflector  was  operated  over  the  frequency  range  24-192  Kcps 
in  the  sea  and  in  the  laboratory.  Excess  attenuations  from  0.15  to 
1.1  db/M  due  to  absorption  and  scattering  by  inhomogeneities  were 
made  over  the  full  frequency  range  with  a  maximum  error  of  less  than 
0.05  db/meter.  Backscattering  at  sea  was  detected  in  the  frequency 
range  24-60  Kcps  with  an  average  scattering  cross-section  per  unit 
volume  of  4.00  cm  /m^  at  a  depth  of  10  feet.  Scattering  was  not 
detected  abbve  60  Kcps.  Variation  of  sound  velocity  with  frequency 
was  measured  from  24-192  Kcps.  The  accuracy  of  these  values  was  0.03 
percent  for  frequencies  over  40  Kcps.  Results  are  compared  and  it  is 
postulated  thatj 

(1)  below  60  Kcps  the  absorption  and  scattering  was  due  to 
microbubbles  of  radius  greater  than  50  microns,  and 

(2)  above  60  Kcps  minute  plankton  caused  absorption  without 
significant  scattering. 
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1.  Introduction 

It  has  been  found  that  the  behavior  of  sound  upon  being 
reflected  from  the  surface  of  the  ocean  does  not  agree  with  theory, 
particularly  for  angles  of  incidence  away  from  the  normal.  One 
of  the  possible  explanations  for  this  disagreement  is  the  presence 
of  small  bubbles  or  other  particulate  matter  near  the  surface. 
Theory  and  laboratory  experience  have  shown  that  a  small  concentra- 
tion of  bubbles  in  water  would  be  particularly  effective  in  causing 
backscattering,  attenuation,  and  changes  in  sound  velocity.  It  is 
the  purpose  of  this  project  to  experimentally  determine  the  presence 
of  microbubbles  and  other  scatterers  and  absorbers  near  the  ocean 
surface . 

A  basic  pulse-echo  system  |_28J  consisting  of  an  open  arrangement 
of  an  electrostatic  piston  transducer  facing  a  reflector  was  used 
to  determine  attenuation,  scattering  and  velocity  dispersion  at  two 
field  locations  described  in  Appendix  VI.  The  three  types  of 
measurements  were  then  used  in  order  to  interpret  the  sources  of 
absorption  and  scattering  of  sound  at  frequencies  ranging  from 
24  Kcps  to  200  Kcps,  in  sea  states  one  and  two,  and  primarily  at 
a  depth  of  ten  feet. 


2.  Instrumentation. 

The  instrumentation  used  in  this  project  consists  of  four 
main  sections;  a  transmitting  section,  a  receiving  section,  a 
transmit-receive  switching  section,  and  a  transducer  section. 
Each  will  be  described  briefly. 

The  transmit  section  takes  a  sinusoidal  signal,  gates  it  to 
form  wave  pulses,  and  then  amplifies  these  wave  pulses.  The  wave 
pulses  are  added  to  a  direct  current  polarizing  voltage  and  then 
applied  to  the  electrostatic  transducer  through  the  transmit- 
receive  switching  section.  Typical  transmitter  outputs  are  shown 
in  figure  2-1.  A  block  diagram  of  the  transmit  section  may  be 
found  in  figure  2-2. 

The  transmit-receive  switching  section  utilizes  a  signal  diode 
network  and  an  electronic  switch  to  isolate  the  transmitter  from 
the  reciever.  Figure  2-3  is  a  block  diagram  of  this  section. 
Figure  2-4  shows  the  receiver  output  while  transmitting  • 

The  receiving  section  consists  of  an  electronic  filter  and 
oscilloscope  display  with  amplification  before  and  after  the  filter. 
Figure  2-5  is  an  example  of  a  receiver  output  display.  Figure  2-^6 
is  a  block  diagram  of  the  receiving  system. 

The  transducer  section  consists  of  a  24  cm  diameter  wide-band 
electrostatic  piston  transducer  and  a  parallel  30  cm  steel 
reflector  located  76  centimeters  along  the  axis  from  the  transducer. 
A  complete  description  of  the  design,  construction  and  characteristics 
of  the  electrostatic  transducer  is  given  in  Appendix  I. 

Figure  2-7  shows  the  electronics  set  up  at  the  Naval  Electronics 


Laboratory  Tower.  Manufacturer's  specifications  for  the  commercial 
equipment  used  in  the  entire  system  may  be  found  in  Appendix  V. 
The  overall  system  characteristics  are  given  in  Appendix  IV, 
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Figure  2-6 
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3 •  Attenuation  Measurements 
a.  Basic  Principle. 

An  Electrostatic  transducer  was  gated  to  emit  a  short 
pulse  and  then  electronically  switched  to  receive  the  successive 
echoes  returned  from  a  reflector  placed  76.4  cms.  in  front  of  it. 
Photographs  of  the  pulse  echo  return  are  shown  in  Fig.  3«1<>  In  all 
cases  the  gain  of  the  oscilloscope  was  adjusted  so  that  the 
amplitude  of  the  second  echo  was  approximately  10  cm.  when  measured 
on  the  oscilloscope  face. 

Successive  echoes  were  reduced  in  amplitude  by  three  types 
of  loss;  divergence,  imperfect  reflection  at  the  reflector  and  the 
transducer  and  attenuation  in  the  medium.  This  attenuation  in  the 
medium  was  due  to  scattering  and  absorption  by  inhomogeneities,  in 
addition  to  the  conventional  absorption  of  the  medium  itself.  The 
object  was  to  measure  the  excess  attenuation  loss  in  the  sea.  For 
all  media  studied,  in  the  frequency  range  24-192  Kcps,  this  medium 
loss  was  smaller  than  the  reference  system  loss  due  to  reflection 
and  divergence. 

At  the  lower  frequencies,  below  approximately  60  Kcps, 
divergence  loss  was  very  great  so  that  the  number  of  pulses  which 
could  be  measured  was  smaller  than  that  desirable  for  graphical  work, 
However  graphs  were  drawn  for  frequencies  as  low  as  24  Kcps  and  used 
for  the  calculation  of  attenuation. 

Fig.  3.1  shows  photographs  from  the  minimum  usable 
frequency  of  24  Kcps  to  the  maximum  frequency  at  which  it  was 
possible  to  operate  the  equipment  satisfactorily,  240  Kcps. 

10 
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a)  24Kcps. 


b) 


AS  Kcps. 


c)  96Kcps.  d)  ,  Kcps> 

.  3.1 
Photographs  ofel  daring  the  calibra  tion  of  the  system. 
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Higher  frequency  operation  will  require  the  use  of  an  amplifier 
capable  of  delivering  an  undistorted  input  to  the  transducer  at 
the  frequency  of  interest  and  a  filter  capable  of  operation  at  that 
frequency.  As  can  be  seen  from  Fig.  3.1  d,  noise  made  the  accurate 
measurement  of  pulse  amplitude  increasingly  difficult  above  200 
Kcps. 

b.  Interference  Effects 

Interference  effects  are  illustrated  by  Fig.  3*2  which 
shows  pulse  amplitude  (logarithmic  scale)  against  echo  number 
(linear  scale).  These  effects  are  superimposed  upon  the  straight 
line  decay  and  increase  in  amplitude  with  echo  number 

In  this  case  the  interference  was  not  due  to  mechanical 
misalignment  which  gives  similar  graphs  but  to  the  development  of 
a  crack  at  the  edge  of  the  transducer.  As  the  aluminum  exposed 
by  the  crack  reacted  with  the  water  and  was  eaten  away,  the 
electrostatically  driven  surface  was  no  longer  circular  and  the 
piston  like  behaviour  was  impaired.  The  interference  effects 
became  worse  with  time  as  can  be  seen  by  comparing  the  1725  and 
1753  curves  beyond  the  tenth  echo. 

c.  Diffraction  Effects  -  Comparison  with  Theory 

Seki,  Granato  and  Truell  [23]  have  used  a  pulse  echo 

system  to  investigate  diffraction  effects  in  the  ultrasonic  field 

of  a  quartz  crystal  acting  as  a  piston  source  and  verify  William's 

[27I  theory.  Our  pulse  echo  system  differs  in  that  the  pulse  is 

repeatedly  reflected  from  a  reflector  the  same  size  as  the  source. 

At  first  thought,  the  method  appears  equivalent  (assuming  perfect 

reflection)  to  the  direct  radiation  from  a  piston.  In  fact,  the 

12 
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physics  more  closely  approachs  the  condition  of  a  sound  beam  that 
is  passed  through  a  series  of  apertures.  At  each  aperture  there 
will  be  a  new  diffraction  (near)  field.  This  problem  has  not 
been  solved,  to  our  knowledge. 

The  comparison  with  William's  theory  is  solely  for  interest. 
It  must  be  emphasised  that  the  determination  of  attenuation  by 
subtracting  the  calibration-slope  from  the  sea-data-slope  bypasses 
the  problem  of  diffraction  that  is  being  discussed  here. 

To  compare  our  field  pattern  with  Seki's  we  have  started 
with  his  data  for  ka  =  50  (curve  a).  This  curve  was  then  modified 
to  include  the  loss  due  to  imperfect  reflection  at  the  transducer 
and  reflector  (Section  I)  and  yielded  curve  b.  It  was  assumed 
that  this  loss  was  independent  of  frequency. 

The  comparison  experimental  data  used  was  that  of  Section 
H,  calibration  of  the  system,  since  at  the  frequency  of  interest 
the  theoretical  absorption  by  fresh  water  was  negligible. 

The  active  radius  of  the  transducer  was  12  cms  so  that, 
at  96  Kcps 

ka  =  2TYf  a 

c 

=  49 

This  is  so  close  to  50  that  the  variation  in  diffraction 
field  due  to  it  is  negligible. 

It  was  also  determined  that  at  96  Kcps  a2/2v_was  93.5  cms 
so  that  the  distance  travelled  by  the  sound  waves  between  echoes 
was  1.64  a.  /  \*     Curve  c  was  then  drawn  with  echo  amplitude  plotted 

14 


against  distance  travelled  by  sound  in  units  of  a2/ \   •  Different 
starting  heights  are  used  for  a)  and  c)  to  preserve  clarity. 

When  d,  the  distance  to  the  first  pulse  of  the  straight 
line  section  of  the  graph  was  plotted  against  a2/X,  for  the 
frequencies  between  24  and  240  Kcps  the  straight  line  of  Fig.  3.4 
was  obtained.  Its  equation  was 

d  -  1.2$  +  13.8  a2/^ 
This  distance  is  the  range  beyond  which  the  attenuation  in  db  is 
proportional  to  the  distance  travelled. 

The  dashed  line  for  2a2/ ^  also  shows  that  diffraction 
effects  and  significant  divergence  losses  occurred  for  beyond  the 
nominal  end  of  the  near  field. 

Diffraction  effects  varied  between  transducers  and  changes 
have  been  observed  in  a  single  transducer  as  shown  by  Fig.  3.5* 
The  additional  noise  in  photograph  3»5a  is  due  to  a  wide  band 
width  setting  on  the  electronic  filter .   Despite  the  changes  in 
early  pulse  amplitudes,  graphs  of  pulse  amplitude  against  echo 
number  yielded  straight  lines  which  were  of  the  same  slope  within 
the  limits  of  accuracy  of  the  experiment  (0.06  db/m). 
d.  System  Alignment. 

To  prevent  interference  effects  of  the  type  shown  in  Fig. 
3.2  it  was  necessary  to  align  the  transducer  and  reflector  so  that 
their  faces  were  as  nearly  parallel  as  possible.  As  shown  in  Fig. 
3.6  the  transducer  mount  could  be  moved  along  the  direction  labelled 
"h"  and  rotated  about  its  supporting  arm,  while  the  reflector  was 
mounted  on  three  adjusting  bolts  so  that  its  direction  might  be 
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adjusted  both  vertically  and  horizontally. 

The  initial  mechanical  alignment  was  made  by  adjusting 
the  transducer  mount  and  the  reflector  until  the  distances  "s"  and 
"d"  were  the  same  between  all  corresponding  points  on  the 
circumferences  of  the  transducer  and  reflector. 

The  system  was  then  placed  in  water  and  operated  at  192 
Kcps  when  gross  misalignment  was  indicated  by  a  decay  pattern  much 

f 

more  irregular  than  those  illustrated  in  Fig.  3.5.  Minor  adjust- 
ments were  then  made  at  this  frequency,  the  most  sensitive  within 
the  normal  operating  range  of  the  equipment. 

When  a  photograph  had  been  taken  and  a  plot  of  pulse 
amplitude  against  frequency  found  to  yield  an  ultimate  straight 
line  as  in  curve  c  of  Fig.  3»3>  the  system  was  regarded  as  being 
satisfactorily  aligned. 

This  procedure  was  difficult  under  laboratory  conditions 
and  proved  more  difficult  while  operating  from  the  N.E.L.  Tower, 
where  the  system  was  mounted  as  shown  in  Fig.  3.7.  Since  the 
transducer  mount  was  easily  rotated,  despite  locking  screws, 
handling  difficulties  prevented  the  maintenance  of  the  accurate 
allignment  needed  for  higher  frequency  operation,  when  working  from 
the  U.S.W.P.G.S.  Oceanographic  Research  Vessel. 

It  is  therefore  recommended  that  the  mounting  should  be 
redesigned  so  that  it  can  be  adjusted  accurately  in  the  laboratory 
and  then  transported  safely  to  the  observation  site. 
e.  Choice  of  Reflector  Material 

To  keep  reflection  loss  as  low  as  possible  reflectors  of 
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steel,  styrofoam  and  corprene  were  made  and  used  in  succession. 
With  a  polarising  voltage  of  200  V.  and  an  input  drive  of  20  V. 
peak  at  60  Kcps  applied  to  the  transducer,  the  amplitudes  of  the 
first  echo  were;   steel,  33  mv.;  styrofoam,  31  mv.j  corprene, 
25  mv. 

The  steel  reflector  was  therefore  selected!  it  had  the 
additional  advantage  of  mechanical  strength  so  that  its  reflection 
coefficient  would  not  be  altered  by  compression  as  its  depth  increased, 

f .  Selection  of  Operating  Voltage 

After  the  reflector  had  been  selected  the  polarising  voltage 
and  input  drive  were  varied  and  the  first  echo  amplitude  observed, 
the  results  being  as  shown  in  Fig,  3. 8.  On  the  basis  of  this  graph 
it  was  decided  to  apply  a  polarising  voltage  of  450  V.,  the  maximum 
which  would  be  accurately  set  with  the  power  supply  in  use,  and 
to  use  an  input  drive  of  up  to  50  V  peak.  This  input  drive  was 
reduced  at  higher  frequencies  to  maintain  a  better  input  pulse  shape. 

g.  Determination  of  Standard  Deviation. 

In  order  to  determine  whether  variations  in  the  field  data 
were  due  to  the  medium  or  the  system  and  to  establish  the  confidence 
level  for  the  calibration,  the  mean  and  standard  deviation  were 
determined  at  60  Kcps  and  192  Kcps.  A  series  of  photographs  was 
taken  as  rapidly  as  possible  to  minimize  changes  in  the  tank  water 
and  equipment.  The  attenuations  were  then  determined  from  each 
photograph. 
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60  Kcps  Result.   (27  samples) 

Mean  X  =  1.62„  db/meter 
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(X1  -  X)2  =  1.58  x  10°2  (db/meter)2 
j=l    J 

Variance  S2  -  6.10  x  10"^  (db/meter)2 

Standard  Deviation  S  =  0.025  db/meter 

192  Kcps  Result. (26  samples) 

Mean  X  =  0.89   db/meter 

Standard  Deviation  S  =  0.015  db/meter 

Fig.  3.9  shows  the  attenuation  in  db/m  vs 

sample  nunber  at  60  kcps.  The  sample  points  are  distributed  in 

mich  the  same  way  as  some  of  the  field  data,  particularly  that  in 

which  depth  and  frequency  were  held  constant  while  photographs  were 

taken  at  2  minute  time  intervals.  The  noise  was  not  significantly 

different  in  the  field.  The  field  curves  have  therefore  been 

smoothed  by  up  to  twice  the  standard  deviation  for  system  noise. 

h.  Calibration  of  System. 

Knowing  the  mean  and  standard  deviation  at  60  and  192 

Kcps,  it  was  desirable  to  determine  the  number  of  photographs  to  be 

taken  at  each  frequency  to  ensure  that  the  error  was  less  than 

1%   at  a  satisfactory  confidence  level.  Since  the  standard  deviation 

was  somewhat  woree  at  192  Kcps  than  at  60  Kcps  all  calculations  were 

performed  at  the  higher  frequency. 

When  n  is  the  number  of  observations 

6  is  the  permissible  error 

s  is  the  standard  deviation 

23 


and  k^  is  from  standard  normal  tables 
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so  that  for  specific  values  of  n,  kji   can  be  determined  and  the 
confidence  level  found. 


NUMBER  OF 
OBSERVATIONS 

n 

CONFIDENCE  1 
LEVEL 

% 

5 
6 

8 

76 
81 
86 

Table  3.1 


From  the  results  shown  in  the  above  table  it  was  decided 
to  take  six  photographs  and  establish  a  confidence  level  of  better 
than  80$. 

In  preparation  for  the  calibration  a  sonar  tank  with  a 
capacity  of  about  10,000  gals,  was  emptied,  washed  down  and  then 
refilled  with  tap  water  to  which  an  algicide  had  been  aded0  After 
filtering  for  one  day  to  remove  particulate  matter,  the  tank  was 
allowed  to  stand  for  four  days  so  that  any  bubbles  of  air  produced 
by  the  filling  might  rise  to  the  surface  or  be  absorbed  in  the  water, 

Six  exposures  were  then  made  at  the  frequencies  of  interest 
and  used  to  establish  the  standard  values  for  attenuation  due  to 
the  system  in  bubble- free  fresh  water. 

At  all  frequencies  the  pulse  amplitudes  for  the  first  and 
last  photographs  were  measured  from  pulse  2  until  they  became  too 
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small  for  valid  measurement  or  the  last  pulse  on  the  photograph 
had  been  measured.  These  pulse  amplitudes  were  then  plotted 
against  echo  number  as  before  to  determine  the  straight  line  section 
of  the  graph  after  which  measurements  from  the  other  four  photo- 
graphs were  confined  to  this  section  of  the  graph. 

From  graphs  of  the  form  shown  in  Fig.  3.3,  curve  c  it  was 
possible  to  deteriaine  the  pulse  amplitude  Pn  at  the  beginning  of  the 
straight  line  or  extrapolated  beyond  it  and  a  second  pulse  amplitude 
Pn+TH  ,  m  bounces  beyond  Pn  but  still  large  enough  to  be  measured 

accurately.  Then  loss  m  bounces  =  20  log    P  n     d-b. 

P  m+n 
From  this  the  loss  in  db/bounce  was  calculated  and  then  converted 

to  db/meter.  knowing  that  the  separation  between  the  transducer  and 

reflector  was  76.4  cms. 

The  results  of  this  determination  yielded  the  curve  shown 

in  Fig.  3.10.  Apparently  algae  were  still  causing  excess  attenuation 

centered  about  1,50  Kcps.  So  that  excess  attenuations  could  be 

determined  from  the  field  data  the  dashed  line  was  drawn  in  and 

has  been  treated  as  the  standard  reference  attenuation  due  to 

divergence  and  reflection  loss. 

i.  Reflection  Losses. 

If  the  attenuation  loss  at  240  Kcps  is  used  to  measure  the 
reflection  coefficient  of  the  system,  divergence  losses  may  be 
neglected  since  the  linear  section  of  the  graph  used  to  calculate 
attenuations  corresponds  to  great  distances  from  the  source. 

Minimum  measured  attenuation  =0.87  db/bounce 

where  the  path  length  for  one  bounce  is  1.53  meters. 

At  the  frequency  of  240  Kcps 
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Absorption  due  to  fresh  water   0.02  db/meter 

Absorption/bounce  due  to  fresh  water  -  0.03  db/meter 

Total  Reflection  Loss  =  0.84  db/bounce 

10  log10  -L.-L   =0.84 
1  2 

System  Reflection  Coefficient  =  -f"""2 

1 

-  0.S3. 
Theoretical  Reflection  Coefficient  of  steel  is  given  by 

Steel  =(/f2^2_^fl_Cl.2 

=  0.94 

steel  X     transducer   =  — - 

II 


since 


Reflection  Coefficient  of  transducer  =  0.88. 

This  determination  has  only  been  made  for  one  transducer. 

If  the  reflection  coefficients  of  other  transducers  differ  from  this, 

erroneous  values  of  excess  attenuation  will  be  found  when  the 

standard  attenuation  is  subtracted  from  the  field  determination  of 

total  attenuation. 

j.  Fresh  Water  Algae  Experiment 

The  anechoic  tank  was  drained,  washed  down  and  refilled  and 

then  allowed  to  stand  for  5  days  without  filtering  or  chemical 

treatment  so  that  the  algae  present  could  develop.  Previous 

absorption  experiments  had  shown  evidence  of  bubbles  after  two  days 

standing  so  that  this  waiting  period  also  allowed  bubbles  to  escape. 

Although  considerable  attenuation  was  measured  between  120 

and  200  Kcps  as  shown  in  Fig.  3.11*  no  back-scattering  (Section  4) 
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was  detected  nor  did  velocity  measurements  (Section  5)  indicate 
any  dispersion. 

A  sample  was  taken  with  a  plankton  net  and  an  analysis 
and  count  kindly  performed  by  Dr.  E.  C.  Haderlie  of  Monterey 
Peninsula  College,  the  algae  population  being: 

Desmids  (20-40  juu)   population  =  20  x  loVliter 

Ceratium  Dinoflagellates  (~90jlo)  population      =  0.4/liter 
Volvox  (100j>u)  population  =  0.2/liter 

This  experiment  showed  that  phytq-plankton  were  capable  of 
producing  considerable  attenuation  without  detectable  scattering  or 
dispersion  at  selected  frequencies.  Whether  there  were  minute 
bubbles  associated  with  the  photosynthesis  of  the  plankton  could 
not  be  determined. 

k.  Absorption  in  San  Diego  Bay 

Attenuation  photographs  were  obtained  at  known  times  and 
depths.  Using  the  Precision  Scientific  Company  Galvanic  Bell 
Oxygen  Analyser  simultaneous  observations  of  temperature  and 
oxygen  content  of  the  water  were  obtained  from  a  thermistor  and 
oxygen  electrode  mounted  on  the  equipment  as  shown  in  Fig,  3.7# 
The  Winkler  tests  needed  to  calibrate  the  oxygen  electrode  were 
kindly  performed  by  Mr.  F,  Crowe  and  Mr.  J.  T.  Hart  of  the  Scripps 
Field  Annex. 

The  photographs  were  analysed  and  graphs  drawn 
to  determine  the  total  attenuation.  From  these  results  the 
attenuations  due  to  the  system  were  subtracted  leaving  the 
attenuation  due  to  the  medium,  as  shown  in  Figs.  3.12  to  3»17« 
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There  was  little  change  in  attenuation  over  a  30  min. 
period  when  observations  were  made  at  a  fixed  frequency. 

When  the  frequency  was  fixed  and  the  depth  varied  a 
consistent  attenuation  pattern  was  produced,  which  appeared  to 
move  slowly  with  time.  Additional  observations  over  24  hr.  periods 
are  needed  to  confirm  that  a  diurnal  change  was  being  observed. 

Small  bubbles  may  be  generated  during  the  daylight  hours 
causing  the  attenuation  at  frequencies  above  96  Kcps  to  increase 
during  the  day.  After  sunset  these  bubbles  would  no  longer  be 
generated  but  would  be  absorbed  into  the  water  or  coalesce  into 
larger  bubbles  so  reducing  high  frequency  attenuation  (See  Fig. 
3.16)  •  Their  coalescence  into  larger  bubbles  would  increase 
attenuation  below  96  Kcps.  as  shown  in  Figs  3»13  and  3»14>  where 
the  evening  attenuations  are  greater  than  the  morning  attenuations. 

At  152  Kcps  and  192  Kcps  there  was  a  definite  relationship 
between  absorption  and  depth,  with  maximum  absorption  near  the 
surface  suggesting  that  the  absorbers  were  organisms  requiring 
light  for  photosynthesis. 

Attenuation  did  not  correlate  with  either  temperature  or 
oxygen  content  of  the  water.  However, the  wetter  was  isothermal 
(in  general  +  0.05°C  or  +  0.10°C)  and  the  oxygen  content  varied 
little  from  100$  saturation  so  that  these  quantities  could  have 
had  little  effect. 

1.  Comparison  of  Attenuation  Values  with  Published  Data. 

In  the  experimental  work  a  short  sound  pulse  passed  many 

times  through  a  small  volume  of  water  so  that  the  measured  attenuation 

was  that  of  a  very  short  path  length  and  was  independent  of  depth, 
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surface  and  bottom  effects.  The  data  published  in  Physics  of  Sound 
in  the  Sea,  Part  1,  Transmission,  measures  loss  over  thousands  of 
yards  where  the  attenuation  was  effected  by  the  surface  and  the 
bottom,  temperature  gradients  and  the  relative  depth  of  the 
transmitter  and  hydrophone  in  use. 

If  Fig.  14,  P102,  graph  of  average  transmission  anomalies 
in  isothermal  water  is  examined  the  attenuation  of  24  Kcps  is 

4.5  db/kyd  =  0.0049db/m  and  the  attenuation  at  60  Kcps  is 
12.5  db/kyd  =  0.014  db/m. 

Similar  values  of  attenuation  are  given  in  Kinsler  and 
Frey,  "Fundamentals  of  Acoustics". 

From  Fig.  3.12  from  data  at  a  depth  of  10  ft.  the  average 
value  of  attenuation  at  24  Kcps  was  1.1  db/m.  about  220  times  the 
accepted  value  and  at  60  Kcps  it  was  0.35  db/m,  only  25  times  the 
accepted  figure. 

Our  data  was  obtained  in  an  environment  which  would  have 
been  unsuitable  for  longer  path  transmission  loss  determinations 
because  of  its  shallowness.  If  our  near  surface  operation  was  made 
in  water  with  many  more  bubbles  and  much  more  plankton  than  would 
be  found  in  a  deep  water  experiment  abnormally  high  attenuations 
would  certainly  result. 

The  only  previous  work  giving  attenuations  of  the  same  order 
of  magnitude  is  that  of  Hansen  and  Barham  \_  9  J  who  used  tanks  resonant 
at  7.9  Kcps  in  the  laboratory  and  2.7  Kcps  in  the  sea  to  determine 
attenuation  by  measuring  reverberation  time.  Working  in  "red  water" 
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areas  they  found  attenuations  greater  than  the  maximum  measurable 
with  their  equipment,  0.405  db/m  at  both  frequencies.  Although 
the  order  of  magnitude  of  their  results  agrees  with  ours  it  must 
be  pointed  out  that 

a)  their  equipment  measures  absorption  not  scattering  and  absorption 

b)  their  equipment  measured  absorption  for  plankton,  not  for 
bubbles  and  plankton. 

Other  possible  explanations  of  our  "high"  value  of 
absorption  have  been  considered? 

1)  Cavitation  loss:  This  is  rejected  because  our  maximum  S.P.L. 
was  70  db  re  1  ubar  whereas  cavitation  would  require  S.P.L.  of 
120  db.  re  1  ubar  at  the  surface  and  greater  than  120  db.  re 

1  ubar  as  depth  increased 0 

2)  There  may  have  been  undetected  bubbles  attached  to  our  trans- 
ducer or  reflector.  This  possibility  cannot  be  completely 
rejected,  although  the  self -consistency  of  our  data  suggests  that 
it  is  an  unlikely  explanation. 
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4o  Detection  of  Baekscatter, 
a0  Theory  0 

A  great  deal  of  work  has  been  done  in  the  study  of  scatterers 
in  a  sound  field.  The  scattering  cross-sections  (  £/<.  )  of  sip, 
solid  or  liquid  particles  is  much  less  than  the  geometrical  cross- 
section  for  kR^  1^ 

kR  =  2TTR/^ 
where    R  =  Radius  of  the  scatterer. 
*  —  wave  length  of  the  soundo 

This  represents  the  well  known  condition  of  Rayleigh 
scattering  where  the  scattered  intesity  is  proportional  to  the 
fourth  power  of  the  frequency. 

The  scattering  cross-section  also  depends  on  the  density 
and  elasticity  of  the  particle  and  the  mediums  the  character!,- 1 
impedance  of  the  scatter  material  (ftC6)  must  be  different  from 
that  of  the  medium  for  scattering  to  occur. 

Figure  4«1  illustrates  the  variation  of  (     h  kR  for  the 
two  extreme  cases  of  heavy j,  rigid  sphere  and  an  ideal  air  bubble. 

When  kR  is  greater  than  2a   both  curves  coincide  and  si- 
oscillations  about  the  mean  falue  of  unity  for  ratic      R  , 
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Figure  4.1 
Comparison  of 
Scattering  and 
Geometrical  Cross- 
Sectionf. 


It  was  assumed  that  only  those  bubbles  at  or  near 
resonance  are  the  cause  of  any  significant  backscatter  at  any 
particular  frequency.  This  is  so  because  the  Q  of  the  bubble  is 
of  the  order  of  10  and  the  effective  scattering  cross-section  at 
resonance  is  enormously  greater  than  the  geometrical  cross-section 
of  the  bubble 
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Where  S  is  the  bubble  damping  constant  due  to  shear  vecosity, 
thermal  conductivity  and  resadiation  as  described  by  Devin  [41  • 
This  relative  scattering  cross-section  is  plotted  in  Figure  4.2  and 
varies  from  918  to  178  at  frequencies  of  20  Kcps  and  200  Kcps 
respectively,  the  frequency  range  of  the  apparatus. 
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Figure  4.6 
Constructed  Backscatter  from  Bubbles. 

Figure  4.5(a)  shows  the  region  with  system  reverberant  noise 
only  while  Figure  4.5(b)  shows  the  increased  amplitude  detected 
when  bubbles  were  introduced.  When  the  noise  is  subtracted  from 
the  signal- plus-  noise  the  difference  is  a  measure  of  the  back- 
scattered  energy  from  the  bubbles;  this  is  plotted  in  Figure  4.6. 

The  photos  of  Figure  4.5  serves  also  to  indicate  the 
difficulty  of  extracting  bubble  scatter  evidence  in  the  presence  of 
reverberant  scatter  from  the  anechoic  tank  walls  as  well  as  the 
system  supports. 

An  experiment  was  conducted  in  the  laboratory  using  a  1.5 
inch  diarnenter  styrofoam  sphere  to  determine  the  effect  of  scatter 
position  on  signal.  The  experiment  showed  that  the  signal  level 
changed  by  no  more  than  1C#  as  the  sphere  was  moved  throughout  the  volume, 
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thus  indicating  the  relative  uniformity  of  the  field.  For  this 
sphere  and  the  frequency  used  the  resulting  ka  of  4  indicated  that 
the  scattering  cross-section  should  be  approximately  equal  to  the 
geometrical  cross-section;  the  experiment  verified  this  and  confirmed 
the  validity  of  the  technique  used. 

The  scatterer  position  experiment  is  described  in  Appendix 
XI. 

c.  Data  Processing. 

1.  Correction  for  Electronic  Noise. 

The  backscatter  is  to  be  evaluated  from  the  amplitude 

H  of  Figure  4.7 •  First,  however,  it  is  necessary  to  use  the  average 

of  the  reflector  echos  A  and  B  (Figure  4.7 )  as  a  measure  of  the 

amplitude  which  would  be  obtained  if  a  "perfect"  reflector  were  at 

the  midpoint  of  the  system. 

Hence : 

A  +  B 

— - —  =  R  «   the  interpolated  echo  level  for  a  perfect 

plane  reflector. 

In  a  similar  manner  the  amplitude  of  H  midway  between  the 

two  echos  can  be  associated  with  the  amplitude  of  the  scatter, 

that  is: 

H-^  Ss.   =  the  total  scattering  cross-section  in  the 

insonified  volume,  including  the  equivalent 
scattering  cross-section  of  the  system  noise. 

Assuming  incoherent  signal  and  noise  we  have; 

H2  =  Eo2  +  N2 

Eo  =  /  H2  -  N2  (4.1) 

46 


wher  1  . 

(3-53  nvj  -  L  voltage. 


.     I   U  at  J/(   ;.: 
jince    -  f  l:. j   r  ■;'.    ct<  r  (o)  j        i 

proportional i 


Sc  it !. 

ithin   ■>..  !-•   .-ir   j~^  ....•- 
R    .     L  .    Area 


c*. 


Ltage\  2 

.-cr.oj 
J 


or 


m 


it     oC 


(f) 


(4.2) 


°st  is  the   bota  .  -section  of  the 

transducer/reflector   tohusae   (V);   to  non   ilj        I     Ls  to  a    mit 
voluiae  °st   L        Lvided   by    /j 


>st 


h7 


Ss.|f^2 


(4.3) 


In  this  investigation  SQ  is  expressed  in  Cm^/m^. 
20  Correction  for  Divergent  Field. 

The  fact  that  the  transducer  radiation  pattern 
changes  with  frequency  results  in  the  insonified  volume  being  a 
variable  as  the  frequency  is  changed. 

Consider  the  situation  as  shown  in  Figure  4.8 


TRANSDUCER 


REFLECTOR 


Figure  4.8.  Divergent  Field  Situation. 
A  pulse  leaving  the  transducer  backscatters  as  it  sweeps 
over  the  truncated  conical  volume  created  by  the  divergence  angle  % 
(for  this  experiment   3°  <  ©  <  9°).  The  fractional  increase 
in  volume  over  that  of  the  original  cylindrical  volume  becomes  s. 


48 


*£**    §d[2rCa+§]_    defane[e7rfa+sffan9j 


—    £ 


Tra^d  paV 


4¥^    rac/af<znQ 


V        a  (4.4) 

The  reflected  pulse  approaching  the  transducer 
also  diverges  and  causes  forward  scatter  from  bubbles;  the 
forward  scatter  from  the  cylindrical  volume  will  be  hidden  in  the 
main  echo  but  the  forward  scatter  from  outside  the  cylinder  will 
arrive  later  due  to  the  longer  path  length  and  will  appear  between 
echos.  This  devergent  region  from  the  reflector  is  subtended  by 
the  same  angle  0  thereby  doubling  the  factional  volume  increases 

41     ~   2d  tan  6  (4.5) 

7      a 


The  ratio  of  effective  folume  to  cylinder  volume 


becomes  i 


l'  "  1  +  %  tan  e  (4.6) 

The  angle  0  is  taken  to  that  of  the  half  power 

half  beam  width  and  is  given  by; 

6  =  sin=1  0.0212 A  ^   0.0212 A  (4.7) 

Because  backscatter  actually  comes  from  the  larger 

volume  V  the  scattering  cross-section  per  unit  volume  requires  a 

correction  factor  of % 
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and  (4.3)  becomes: 


Ss  = 


_  S 


Kx  Eo 


^   R  j  (4.8) 

The  range  of  Divergent  Field  connection  factors  is 
given  in  Table  4.1  for  frequencies  between  20  Kcps  and  60  Kcps. 
Table  4.1  Divergent  Field  Connection  Factors  (Kj_) 


Frequency  (Kcps) 

20 
30 
40 
50 
60 


Kl 

0.600 
O.676 
0.730 
0.764 
0.794 


3.  Correction  for  Pulse  Duration. 

Consider  the  two  successive  pulse  as  in  Figure  4.9 


2°=  0.5 msec 
7~=  /.oa  msec 


Figure  4.9  Pulse  Correction  Situation 
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During  the  time  of  the  pulse  (f )  any  backscattered 
energy  will  be  obscured  by  the  much  greater  level  of  the  pulse, 
to  correct  for  this  a  comparison  is  made  between  energy  levels 
rather  than  power  levels,  assuming  that  the  scatter  power  is 
spread  over  the  full  time  between  the  pulses  (T.) 

Then  4.8  becomes: 


In  general  for  this  experiment i 


-£- 


so     S.  =      S     fj3j2jo 


's 


f   Ki  K2  Ep    \ 


2 


(4.10) 
4.  Correction  for  Isotropic  Scattering 

Because  the  scattering  from  any  resonant  bubbles 
will  be  isotropic  in  nature  the  transducer  will  only  intercept  a 
portion  of  the  spherical  scatter  in  an  amount  depending  on  the 
position  of  the  bubble  within  the  transducer/reflector  volume. 

In  order  to  obtain  a  correction  factor  for  this 
effect  it  was  necessary  to  obtain  a  value  for  the  mean  solid  angle 
subtended  at  the  transducer  assuming  a  homogeneous  distribution 
of  bubbles  within  the  volume. 

Consider  a  bubble  at  position  P  in  Figure  4.10. 
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Figure  4.10 
Geometry  of  Isotropic  Scattering  Correction  Factor 

The  solid  angle  XI  is  given  by: 

H.  -  2TT  (1-cos^)  (4.11) 

and 

Cos( 


4  «4  +  £osJi 


From  the  Law  of  Cosines: 


Cos  2+  =  a2  +  b2  -  °2 
T        2ab 


where:  a2  =  (R-r)2  +  X2 

c2  =  4R2 

b2  =  X2  +  (R+r)2 
hence:  XI  =  2  TT  (l  -   /  \  +   Cos  2(j)"  )     (4.12) 


expanding : 
Cos2(j)  = 


X2  +  r2  -  R2 


^1^X2  -  27.2rX  +  X4  +  342OO  -  370r2  +  r^  + 
..52 


-■-i  — (4-L3 ) 

+27.2  rX2  +  r2X^  +  3^0  X2 


Cos  2$  =  2B 

A 


-*r*-fi       A 


fl  =2-tt  (1  -  /*+  -JL  1  (4.14) 


To  find  the  mean  solid  angle  fl   it  is  necessary 
to  perform  a  double  integration  in  the  two  variables  r  and  X: 

R    L 

ft  =  2TT  I         /    -flLrdr  dx 
J      J         TT  R2L 
r=0     X=0 

R         L 

R2l    7     7 

R=0     X=0 

Substituting  the  expression  for  XI  : 
L   R 

X=0  r=0 
Substituting  the  values:  R  =  13.6  cm 

L  =  76.4  cm 
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76.4    13.6 

•fl  -     /i  + 

/I 


*=lk    j    f     '(1-/*7rfO**< 


1=0 
76.4    13.6 


n  =  2ir-nkJ  7    7i  +  7T  drdx    (4-16) 

0    0 
An  approximate  value  was  obtained  for  il  by 
finding  the  mean  value  of  fl   along  the  axis,  (for  the  case  r=0). 
Hence:  A  %  X^  +  34200  +  370X2 

B  «s  (X2  r2)  J 
Substituting  appropriate  values  gave: 

n  %  0.5  TT 
To  find  the  total  scattering  ross-section  per 
unit  volume  it  was  necessary  to  multiply  the  experimentally 
obtained  values  by  Kq  where 

Ko  =  Total  Spherical  Solid  Angle  =  477-  =  g 

■ft  0.5TT 

5.  Tank  Calibration  for  Reverberation. 

For  various  data  rums  Ss  can  be  calculated  and 
plotted  as  a  function  of  frequency.  In  general,  however,  there  was 
an  unknown  amount  of  interference  due  to  the  transducer  support 
structure.  An  attempt  to  correct  for  this  a  number  of  calibration 
runs  were  made  in  the  USNPGS  ancohoic  tank.  It  was  found  that  most 
of  the  reverberent  level  measured  in  the  anechoic  tank  was  due  to 
scattering  from  the  tank  itself.  Orientation  of  the  transducer- 
reflector  had  a  substantial  effect  on  background  and  so  the  curve 
of  background  reverberation  obtained  in  the  "calibration"  over- 
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corrects.  The  magnitude  of  the  correction  curve  was  considerably 
less  than  the  signal  levels  and  any  correction  was  less  than  10$. 
6.  Summary  of  Corrections  for  Backscatter  Analysis. 

Combining  all  the  correction  factors  (4.10)  becomes: 

Ss  -  -4*3  (Jf*-) 
finally:  SS=-^(J^L_^         (4.18) 


where     K  =  K3K2  V  K3 

Values  for  K  over  the  frequency  range  20-60  Kcps 
are  given  in  Table  4.2. 

T^ble  4.2  Total  Correction  Factor  (K). 
Frequency  (Kcps)  K 

20  2.42 

30  2.73 

40  2.95 

50  3.09 

60  3.21 

It  must  be  conceeded  that  the  number  and  magnitude 

of  the  corrections  that  have  been  discussed  above  lend  uncertainty 

to  the  final  calculation  of  scattering  cross-section  per  unit  volume. 

The  writer  feels  that  the  final  correction  itself  may  be  in  error 

by  a  factor  perhaps  as  great  as  two.  This  error  factor  will 

reflect  in  all  values  of  scattering  cross-sections  which  follow. 

d.  Analysis  of  Data 

An  analyis  is  now  carried  out  assuming  that  the  scattering 

cross-sections  obtained  were  due  only  to  the  presence  of  bubbles; 

this  permits  calculation  of  the  maximum  bubble  density  and  volume 

concentration  for  any  particular  datum. 
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corrects.  The  magnitude  of  the  correction  curve  was  considerably 
less  than  the  signal  levels  and  any  correction  was  less  than  10$. 
6,  Summary  of  Corrections  for  Backscatter  Analysis. 

Combining  all  the  correction  factors  (4.10)  becomes: 
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where  K  =  K3K2 {   E3 

Values  for  K  over  the  frequency  range  20-60  Kcps 
are  given  in  Table  4.2. 

Tqble  4.2  Total  Correction  Factor  (K). 
Frequency  (Kcps)  K 

20  2.42 

30  2.73 

40  2.95 

50  3.09 

60  3.21 

It  must  be  conceeded  that  the  number  and  magnitude 

of  the  corrections  that  have  been  discussed  above  lend  uncertainty 

to  the  final  calculation  of  scattering  cross-section  per  unit  volume. 

The  writer  feels  that  the  final  correction  itself  may  be  in  error 

by  a  factor  perhaps  as  great  as  two.  This  error  factor  will 

reflect  in  all  values  of  scattering  cross-sections  which  follow. 

d.  Analysis  of  Data 

An  analyis  is  now  carried  out  assuming  that  the  scattering 

cross-sections  obtained  were  due  only  to  the  presence  of  bubbles,; 

this  permits  calculation  of  the  maximum  bubble  density  and  volume 

concentration  for  any  particular  datum. 
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M. 


From  classical  theory  I  22  J  ,  the  total  scattering  cross- 
section  per  unit  volume  is  given  by: 

oo 


4  f    °WJ>  (4.19) 


s 


/W-1)** 


Where  n(R)dR  is  the  number  of  bubbles  per  unit  volume 
which  have  a  radius  between  R  and  R+dr. 

An  approximate  evaluation  of  S3  is  carried  out  by 
substituting  the  small  quantity. 

Q    ^Wo_i  =  _R_-landdft=_dR 
^   Wo       Ro  l   Ro 

and  integrating; 


(4.20) 


The  integrand  is  simplified  by  dropping  higher  order  terms 
of  the  small  quantity  R  so  that  we  need  only  consider: 

R 

dR *JLtan-l  -^ -A       (4.21) 


-k 


if  we  assui  ie  •>    ss  0.1 
then  A  «,  10  tan"1  20R 
Amax  =  15.7 
if  we  let  R  =  0.05 
A^7.S3 


and  __A_  =  LJL.  «-   0.50 
Amax    15.7 
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In  other  words,  the  bubbles  within  5%   of  the  resonant 
size  will  cause  50$  of  the  observed  scattering  cross-section  per 
unit  volume  (Ss). 

The  5%   bandwidth  was  used  in  the  determination  of  the  number 
of  bubbles  per  unit  volume  with  radii  between  r  and  r  +  dr  (hereafter 
called  "bubble  density"  =  n(R)dR)  and  the  ratio  of  bubble  air 
volume  to  water  valume  (hereafter  called  "bubble  volume 
concentration"  =  U(R)dR). 

The  bubble  density  over  any  of  the  5%   increments  of 
frequency  is: 


n(r)dR  =      2  (4.22) 

S   (Wo) 

where  ©^(WO)  is  the  resonant  scattering  cross-section  of  a 
single  bubble  as  shown  in  Figure  4.3. 


6>S(W0)  = 

4-n-Ro* 
S2 

Hence 

n(R)dr  =  - 

*s  S,2 

(4.23) 

The  bubble  volume  concentration  over  the  same  frequency 
increment  is: 

U(R)dr  =  -4-TT  Ro3  n(R)dr         (4.24) 

Data  was  taken  in  San  Diego  Bay  using  the  USN-NEL 
oceanographic  Research  Tower  as  an  equipment  platform  and  in 
Monterey  Bay  using  the  USNPGS  Hydrographic  Research  Vessel. 
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The  variables  considered  in  gathering  the  scattering  data 
were:  depth,  frequency,  time,  day /night  comparison  and  sea  state. 
It  was  found  that  the  system  was  frequency  limited  to  a  range  of 
20  Kcps  to  60  Kcps  when  using  this  technique;  the  oscillator  had 
a  lower  limit  of  ~~  20  Kcps  and  above  60  Kcps  the  signal  to  noise 
ration  deteriorated  to  a  point  of  unusability. 

Most  of  the  data  were  gathered  at  a  depth  of  10  feet 
which  seemed  to  be  the  optimum  depth  for  scatter  detection  without 
being  influenced  by  wave  action  and  surface  effects. 

In  measuring  the  signal  level  it  was  found  that  there  were 
interference  or  reverberation  effects  associated  with  the 
configuration  of  the  apparatus;  this  is  discussed  in  Section  4»E. 

The  interference  resulted  in  a  standing  wave  type  pattern 
being  superimposed  on  the  signal  due  to  scattering;  becuase  of  this 
the  measuring  technique  used  was  to  find  and  measure  all  the  signal 
mayi ma  and  minima  as  frequency  was  varied;  when  all  the  points 
were  plotted  against  frequency  ar  smoothing  technique  was  used  to 
approximate  the  signal  due  to  scattering  alone.  Examples  of  the 
field  data  and  related  "smooth"  curve  are  shown  in  Figure  4.11. 

The  "smooth"  curve  is  used  as  the  basis  for  subsequent 
curves . 

e.  Interference  Problems. 

All  the  data  gathered  from  the  NEL  Oceanographic  Research 
Tower  suffered  from  spurious  scattered  radiation.  Leaks  of 
interference  were  spaced  at  about  5  Kcps  intervals;  this  indicated 
the  possibility  of  an  interfering  source  about  1  foot  from  the 
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latter  observation  might  be  explainable  by  phytoplankton 
photosynthesis  during  the  day  with  emission  of  the  smaller  bubbles 
which  coalesce  after  dark,  decreasing  the  number  of  small  bubbles 
and  increasing  the  nunber  of  large  bubbles  which  rise  during  the 
night  and  disappear  with  daylight  again;  there  might  also  be  a 
vertical  movement  of  bubbles  fixed  to  plankton  during  the  night. 
A  24  hour  observation  at  40Kcps  (*igure  4.21)  with  readings 
every  hour  confirmed  the  increased  scattering  activity  during  the 
dark  hours  at  this  frequency;  the  scattering  level  rose  sharply 
between  2000  and  0800  hours. 

Assuming  that  the  scattering  is  due  solely  to  bubbles  the 
"bubble"  density  curves  for  runs  A  to  E  are  plotted  in  Figure  4.22 
to  4.24.  In  all  cases  the  density  distributions  peak  between  38 
Kcps  and  47  Kcps;  this  is  in  general  agreement  with  the  work  done 
by  Glotov,  Kolobaev  and  Neuimim  I  6 J  in  which  it  was  observed  that 
the  distribution  of  wind-initiated  bubbles  peaked  at  45  Kcps.  The 
comparison  may  even  be  better  when  one  considers  that  the  Glotov 
experiment  was  done  in  a  tank  with  "clean"  sea  water  whereas  the 
investigation  of  this  paper  conserns  bubbles  in  the  ocean  which 
may  be  attached  to  particulate  matter  which  would  tend  to  increase 
the  bubble  mass,  reducing  the  resonant  frequency  and  shifting  the 
distribution  peak  to  a  lower  frequency. 

There  is  approximately  a  4 si  ratio  between  the  highest  and 
lowest  density  peaks  observed. 

The  corresponding  "bubble"  volume  concentration  for  runs 
A  to  E  are  shown  in  Figures  4.25  to  4.27. 
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A  number  of  runs  (F  to  L)  were  conducted  to  determine  the 
variation  of  back-scatter  with  depth.  Figures  4.28  and  4.29  show 
comparitive  day /night  data  at  40  Kcps;  in  general  there  is  a  drop- 
off with  increasing  depth.  The  drop-off  at  night  is  less  pronounced 
and  the  night  scatter  is  less  near  the  surface.  At  50  Kcps  however 
(Figure  4.30)  the  night  scatter  level  is  greater  than  the  day  level 
at  all  depths  and  in  both  cases  the  level  is  essentially  unchanged 
with  depth  indicating  that  bubbles  of  this  size  (radius  =64  microns) 
are  more  evenly  distributed  in  the  medium  then  the  40  Kcps  bubbles 
(radius  =  80  microns)  which  are  more  concentrated  near  the  surface. 

The  depth  run  in  Monterey  Bay  (Figure  4.31)  indicated  that 
there  was  little  change  in  the  scattering  level  to  a  depth  of  56 
feet  at  a  frequency  of  35.7  Kcps.  However,  there  was  no  data 
between  the  surface  and  10  feet  where  the  level  was  probably  higher. 
Another  reason  for  the  apparent  lack  of  scattering  drop  off  with 
depth  is  the  fact  that  these  data  were  taken  in  a  condition  of  sea 
state  1  whereas  most  of  the  previous  data  were  for  sea  state  2 
conditions  in  which  it  would  be  expected  that  there  would  be  more 
activity  near  the  surface,  particularly  at  such  a  low  frequency. 

Referring  to  the  total  scattering  cross-sections  per  unit 
volume  as  plotted  in  Figures  4.18  to  4.20  the  mean  attenuation  due 
to  the  "bubbles -only"  assumption  is  approximately  1.0  db  per  meter 
in  the  frequency  range  30  to  60  Kcps;  this  figure  is  obtained  by 
using  the  expression  jj22j  ; 

a 


=  4.34/A.jSs   db/meter 
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5.  Velocity  Measurements, 
a.  Brief  Theory. 

The  basi  •  properties  of  a  liquid  that  control  the  velocity 
of  sound  through  the  liquid  are  compressibility ,  ^  and  density,,  p. 
The  sound  velocity  can  be  expressed  as 

c  =  1/  /X^  (5-1) 

This  equation  can  be  derived  from  the  equation  of  state  f*21  J  . 

The  plane  wave  solution  to  the  one-dimensional  wave  equation s 


can  be  written 


|  =  Ae"^  eJw* 
where  '         k  =  c*  +  j/8  =  §L  (5-2) 

It  can  be  seen  that  o<  will  describe  the  attenuation  of  the  wave  with 
increasing  xs   and  0,  will  describe  the  propagation  constant  of  the 
wave. 

If  we  substitute  a  complex  Xj> 
K  =  a  +  jb 
into  equation  5=1  and  then  combine  equation  5-1  with  5-2  we  obtain 

c  =  ^7  <«> 

For  bubble  free  water 


a 


VyP-V 


where  c0  is  the  sound  velocity  in  bubble  free  water. 

Since  we  are  looking  for  the  effect  of  small  bubbles  in  small 
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concentrations  we  will  assume  that  the  density  of  water  with 
bubbles  is  not  different  from  that  without  bubbles.  With  this 
in  mind,  if  we  have  an  equation  for  the  compressibility  of  water 
with  bubbles,  we  can  calculate  the  velocity  of  sound  through  bubbly 
water.  Meyer  and  Skudrzyk  have  derived  such  an  equation  for 
bubbles  of  a  particular  size  L20J  • 


R?%2 


2 
o 


< 
to 


-i-  -  i  +  j  d  JL 


where 

V  =  The  volume  concentration  of  bubbles  in  water  =  ratio  of 

volume  of  air  in  bubbles  to  volume  of  water 
O  =  water  density 
1^0     =  bubble  resonant  frequency 

£0  =  frequency  (independant  variable) 
R  =  bubble  radius 
d  =  bubble  damping  factor 

)\,  =  added  compressibility  due  to  bubbles  in  water. 
The  real  part  of  )C  is 


a7  --2L 


/O^Ufc- 


%' 


-  1 


f--; 


Ud2  **>* 


«u* 


(5-4) 
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The  velocity  in  bubble  water  is  then 


co 


c  = 


/(a+a')/0 


%2 


-~y—  -   1 
1  +  ^JI.  -J^L 1 


(£->)s* 


"XJ   (5-5) 

Fox  Curley,  and  Larsen  have  obtained  experimental  data  that  qualita- 
tively confirm  Meyer  and  Skudrzyk's  equation  for  a  distribution 
of  bubbles  centered  at  60  microns  radius  and  ranging  from  20  to 
120  microns  radius.  This  bubble  distribution  constituted  0.02 
percent  of  the  water  volume,  |_6  J 

If  we  know  the  volume  distribution  of  bubbles  over  a  range 
of  sizes  we  can  integrate  equation  5-4  over  R  (note  that  ^,  ■•'*  (R)) 
and  substitute  the  new  a1  total  into  equation  5-5  to  determine 
velocity  as  a  function  of  frequency. 

This  was  done  using  a  CDC  1604  computer  and  the  results 
for  various  distributions  are  shown  in  figure  5-1 •  The  computer 
program  is  presented  in  Appendix  Vii.  The  volume  distributions  put 
into  the  computer  were  divided  such  that  volume  concentration  values 
were  assigned  to  discrete  bands  of  bubble  sizes.  The  band  width 
used  was  to  microns.  An  alternative  method  of  incrementing  the 
bubble  wolume  distribution  would  be  to  assign  volume  concentration 
values  to  discrete  bands  of  bubble  resonant  frequencies.  The  former 
method  was  used  because  it  appears  more  often  in  the  literature. 
[}>   22,  26] 
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5.2  Experimental  Methods. 

In  order  to  determine  dispersion  due  to  bubbles  in  the  ocean 
the  velocimeter  used  must  meet  the  following  requirements.  First, 
it  must  have  continuous  frequency  coverage  from  10  Kcps  to  200  Kcps; 
these  frequencies  are  determined  from  the  characteristics  of  dis- 
persion curves  shown  in  figure  5-1  that  are  computed  from  the  bubble 
density  distribution  observed  by  Glotov,  Kolobaev  and  Neuimin.  [J>  J 
Second,  it  must  be  quickly  and  easily  read  in  order  to  sweep  the 
frequency  range  before  ocean  conditions  change. 

Any  dispersion  experiment  will  be  limited  by  temperature 
inhomogeneities  in  the  ocean.  Liebermann  has  measured  these 
inhomogeneities  at  depths  of  100-200  feet  and  found  a  typical 
temperature  variation  of  ,05°C  occurring  over  a  correlation 
distance  of  60  centimeters.  This  is  equivalent  to  a  velocity 
variation  of  0.2  M/S  or  0.01  percent  occurring  during  an  average 
period  of  0.6  seconds  assuming  a  current  of  two  knots  flows  by  the 
transducer.  Similar  data  taken  at  the  sea  surface  as  a  function 
of  time  is  not  available.  Figure  5-1  shows  that  if  the  velocity 
error  due  to  temperature  is  of  order  0.2  M/S  the  smallest  total 
volume  concentration  of  bubbles  in  the  Glotov  distribution  that 

could  be  detected  in  ocean  velocity  measurement s  would  be  about 

—8 
6  x  10  -.  .  Salinity  inhomogeneities  could  account  for  a  variation 

of  1.23  M/S  per  part' per  thousand  change  in  salinity  but  whether 

such  changes  can  be  expected  during  an  experiment  is  not  known. 

Of  the  various  methods  for  measuring  sound  velocity  the  time 

of  flight  method  appears  best  suited  for  this  project  in  that 
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frequency  can  be  changed  continuously  without  changing  the  physical 
dimensions  of  the  transducer  assembly  and  the  same  basic  pulse 
echo  system  can  be  used  in  making  attenuation  and  backscattering 
measurements  as  well.  A  variation  on  the  time  of  flight  technique, 
the  "sing  around"  technique  could  not  be  used  because  it  would 
require  major  changes  in  electronic  systems  settings  thereby  pre- 
cluding simultaneous  velocity,  attenuation,  and  backscattering 
measurements.  |_8j   Also  the  "sing  around"  technique  is  not  accurate 
in  a  dispersive  medium.  Time  of  flight  measurements  have  been  made 
in  the  laboratory  by  Greenspan  and  Tschiegg  with  accuracies  of  one 
part  in  30,000  or  .003  percent.  [7] 

Two  different  techniques  for  measuring  time  of  flight  have, 
been  tried.  Each  will  be  described  below. 
Envelope  Detector  and  Echo  Counter 

In  this  technique  the  time  interval  is  measured  from  the  start 
of  the  main  bang  to  the  start  of  a  preselected  echo.  Figure  5-2  is 
a  block  diagram  of  the  system. 

The  advantage  of  this  technique  lies  in  its  ease  of  operation 
and  reading.  The  one  megacycle  counter  gives  an  instantaneous 
reading  of  travel  time.  In  the  laboratory,  travel  time  readings 
to  the  sixth  echo  were  obtained.  This  provided  a  reading  resolution 
of  one  microsecond  in  about  6300  microseconds  or  about  0.016  percent, 
However,  the  reading  was  sensitive  to  frequency  and  echo  amplitude 
changes  as  well  as  velocity  so  that  the  accuracy  of  the  velocity 
value  was  about  5%,     A  discussion  of  this  error  and  a  possible 
remedy  may  be  found  in  Appendix  VIII. 
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Figure  5=2 
Block  diagram  of  envelope  detector  and  echo  counter. 

Time  Marker. 

This  system  provides  an  oscilloscope  display  similar  to  that 
of  a  radar  range  marker  on  an  expanded  scale  "A"  scope.  Figure  5=3 
is  a  block  diagram  of  the  system.  Figure  5-4  is  a  typical  display 
provided  by  this  system.  Appendix  IX  presents  a  detailed  description 
of  the  system. 

The  time  marker  system  appeared  very  attractive  because  of  the 
precision  with  which  measurements  could  be  made.  When  measuring 
travel  time  to  the  11th  eeho  a  reading  resolution  of  one  part  in 
10,000  was  possible. 
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The  method  of  measurement  was  as  follows.  Time  was  measured 
to  an  easily  identifiable  point  on  the  first  echo.  The  time  was 
then  measured  to  the  corresponding  point  on  a  later  echo.  The 
velocity  was  then  obtained  from 

c  —   distance  traveled  from  first  echo  to  later  echo 
travel  time  from  first  echo  to  later  echo 

Difficulty  was  encountered  in  identifying  the  proper  point  on 
the  later  echo  that  corresponded  to  the  previously  selected  point 
on  the  first  echo.  It  is  believed  that  this  difficulty  may  have 
been  caused  by  phase  shift  upon  reflection  from  the  transducer  and 
misalignment  of  the  transducer  reflector  apparatus 0  The  error 
due  to  phase  shift  upon  reflection  was  determined  from  travel  time 
data  taken  in  an  anechoic  tank  filled  with  water  believed  to  be 
non-dispersive^-.  Figure  5=5  is  plotted  from  such  data  and  used  as 
calibration  curves  for  data  taken  using  the  same  transducers. 
It  should  be  noted  that  the  rapid  apparant  change  of  travel  time 
with  frequency  at  low  frequencies  will  cause  low  frequency  determ- 
ination of  travel  time  to  be  inaccurate.  This  effect s   due  to  trans- 
ition from  stiffness  to  resistance  control  in  the  transducer  is 
reflected  as  a  phase  shift  which  is  a  function  of  frequency. 2  Any 
lack  of  parallelism  causes  extra  cycles  to  appear  in  front  of  late 
echoes  at  high  frequencies.  This  ambiguity  in  echo  leading  edge  due 
to  misalignment  is  a  function  of  frequency  which  makes  it  difficult 
to  follow  the  preselected  point  on  late  echoes  as  frequency  is 
See  Section  3H  for  a  discussion  of  tank  preparation . 
^See  Appendix  I  for  a  discussion  of  this  transition. 
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If  the  measuring  point  on  the  later  echo  does  not  correspond 
to  the  measuring  point  on  the  first  echo,  the  measured  velocity- 
plotted  as  function  of  frequency  in  a  non-dispersive  medium  will 
be  of  the  form 


0=  5-rxr  (5"6) 


where 

c  =  Measured  velocity 

d  =  distance  traveled 

To  =  correct  travel  time 

K  =  Error  in  selecting  corresponding  point  on  later  echo 
(in  cycles  or  franctions  of  cycles) 

f  =  signal  frequency 
By  fitting  curves  generated  from  equation  5-6  for  various  values 
of  K  to  the  experimental  curve,  an  approximate  value  for  K  can  be 
obtained  from  the  curve  of  best  fit.  The  measuring  point  on  the 
later  echo  can  now  be  corrected  by  the  amount  K.  Figure  5-6  ill- 
ustrates the  above  discussion.  Lack  of  parallelism  is  not  a  source 
of  error  dice  K  is  determined.  Figure  5-7  shows  how  the  signals 
received  from  all  points  on  the  reflector  will  always  add  to  present 
a  signal  that  represents  an  echo  from  an  effective  parallel  reflector 
located  some  average  distance  from  the  transducer  that  does  not 
change  with  frequency.  The  error  due  to  uncertainty  in  this  distance 
will  cause  an  error  of  0.96  M/S  per  millimeter;  this  error  does  not 
change  with  frequency. 
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Reading  error  may  be  caused  by  instability  in  the  delay  and 
pulse  width  in  the  pulse  generator 5  instability  in  the  oscilloscope 
synchronizing  circuity  instability  in  the  trigger  points  in  the 
counter  and  rise  time  of  the  trailing  edge  of  the  delayed  pulse 
which  was  used  as  the  time  marker.  The  rise  time  error  was 
eliminated  by  consistantly  taking  measurements  with  the  desired 
point  on  the  echo  one  centimeter  below  the  top  of  the  time  marker 
pulse.  Errors  due  to  instability  were  determined  by  taking  many 
readings  at  a  few  selected  points  and  plotting  distribution  curves. 
The  distribution  curves  are  shown  in  figure  IX-1.  The  reading 
error  as  determined  from  these  curves  is  conservatively  estimated 
as  .03  percent o 

The  field  measurements  of  velocity  are  discussed  in  the  follow- 
ing section,. 
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6.  Data  Evaluation. 

A,  a)  Comparison  of  attenuation  and  back-scattering  data* 

To  compare  the  data  obtained  in  Sections  3  and  4  it 
is  convenient  to  convert  both  to  the  equivalent  extinction  cross- 
section  using  the  relationships  developed  in  Physics  of  Sound  in 
the  Sea  [l8]  . 

The  basic  attenuation  relationship  is 
a  =  +  4.34  cre  N(R> 
=  +  4.34  Se, 
where  a  is  the  attenuation  measured  in  db/meter 

(re   is  the  extinction  cross-section  of  a  single  bubble 
N(R)  is  the  number  of  bubbles  of  required  radius  and 
Se  is  the  total  extinction  cross-section. 
Then  using  the  experimental  data 
Se  =  0.23a  m-1  . 
To  determine  Se,  given  Ss,  the  basic  relationship  for 
a  single  bubble  is 

°s  =   KB 

o-e  S 

but  KR  =  0rso  that  replacing  the  single  bubble  case  by  that  for 
multiple  bubbles, 

b 
s°=  X   s" 

where  0  is  the  total  damping  constant 

and   x  is  the  radiation  of  damping  constant. 
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A,  b)  Experimental  Results 

Fig.  6.1a  shows  the  extinction  cross-sections 
calculated  from  scattering  cross-section  and  data  obtained  directly 
from  attenuation  measurement.  In  view  of  the  small  initial  signal 
return  from  the  scatterers  and  the  large  number  of  correction 
factors  applied,  agreement  between  the  curves  is  satisfactory. 

Since  no  scattering  is  observed  above  60  Kcps  it  is 
believed  that  the  effect  of  microbubbles  is  negligible  above  that 
frequency.  This  implies  that  bubbles  of  radius  less  than  50  microns 
are  rapidly  dissolved  while  bubbles  of  greater  than  $0  micron  radius, 
possibly  attached  to  particles,  are  persistent. 

Attenuation  at  frequencies  higher  than  60  Kcps  is 
probably  due  to  some  mechanism  other  than  scattering  and  absorption 
by  microbubbles.  This  is  believed  to  be  absorption  by  minute 
plankton  with  negligible  scattering. 

Fig.  6.1b  shows  extinction  cross-section  as  a  function 
of  frequency  and  indicates  the  frequency  ranges  over  which  the 
effects  of  microbubbles  and  plankton  are  dominant. 
A.  c)  Conditions  for  Absorption  without  Scattering. 

Data  both  from  the  sea  and  from  the  fresh  water  algae 
experiment  have  shown  considerable  attenuation  where  no  scattering 
was  detected^  so  that  bubbles  cannot  play  a  major  part  in  this 
process  at  these  frequencies. 

Meister  and  St.  Laurent  £l9j  have  made  a  study  of  the 
absorption  due  to  algae  in  the  frequency  range  15-27  mcps  and  con- 
cluded that  scattering  was  not  a  mechanism  for  absorption.  They 

also  coj  f 3 u ;j  structural  viscosity       *  he 
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algae  behaved  in  the  same  manner  as  the  structural  viscosity- 
present  in  the  liquid „ 

Lax  and  Feshback  [l3]   show  that  for  particles  suffic- 
iently small  with  respect  to  a  wavelength  absorption  will  be 
dominant  and  scattering  will  be  negligible.   This  occurs  when 
ka<  Ool,  where  a  is  the  particle  radius 

and  k  =  2  7Y  f 
c 

The  maximum  radius  for  absorption  without  scattering  is 


0.1  c 


2  Tf  f 
=  2.4  x  1CA/F  microns 
where  F  is  frequency  in  Kcps. 

This  gives  the  maximum  particle  radius  for  absorption 
with  negligible  scattering  at  the  frequencies  shown  in  the  table. 


Frequency 

a-  max 

Kcps 

microns 

24 

1,000 

30 

800 

38 

630 

48 

500 

60 

400 

76 

320 

96 

250 

120 

200 

152 

160 

192 

125 

Table  6.1 


Both  these  studies  confirm  that  plankton  could  give 

attenuation  without  scattering  over  the  frequency  range  of  the 

measurements 9   but  neither  suggests  a  mechanism  which  would  make 

the  absorption  a  function  of  frequency^  although  definite 
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attenuation  peaks  were  obtained  in  the  fresh  water  algae  experiment, 
B.  a)  Bubble  distribution  from  data. 

To  determine  the  volume  of  air  per  volume  of  water, 
assuming  the  attenuation  to  be  due  to  bubbles,  use  was  again  made 
of  the  N.D.R.C.  relationships  (18  J  • 

If  R,  is  given  by  d  R 

where  RQ  is  the  radius  of  the  bubble  resonant  at  the  frequency  of 
interest,  the  fraction,  p,  of  the  total  extinction  cross-section 
which  is  due  to  bubbles  between  R  -  dr  and  R  +  dr,  is 


1 

tan  -1  2£ 
ft 

n 

28 

= 

2 
Tf 

tan 

2A, 
b 

Then  the  attenuation  due  to  bubbles  in  this  range  is  simply  pa. 

Defining  u(R)  dR  as  the  volume  of  air  per  unit  volume 
for  the  bubbles  under  consideration 

pa  =  1.50  x  10+3  u(r)  db/meter 
u(R)  =  6.67  x  1(H*  pa 
u(R)  dR  =  2  u(r)  dr 

=  1.33  x  10-3  pa  dr. 
Bubble  distribution  and  volume  concentration  bases  on 
scattering  data  were  illustrated  in  detail  in  section  4. f»  but 
certain  curves  are  repeated  here  for  comparison  with  the  distribution 
and  concentrations  obtained  from  the  attenuation  data. 
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Figure  6-.^  shows  bubble  volume  concentration  obtained 
from  scattering  data  in  3an  Diego  Bay  while  figure  6-3  shows  the 
bubble  volume  concentrations  needed  to  produce  the  attenuations 
shown  in  ^igure  3-12 „  The  bubbles  of  radius  less  than  50  microns 
are  not  belived  to  exist  in  significant  amounts. 

Figure  6=4  shows  the  bubble  volume  concentrations 
obtained  from  scattering  data  in  Monterey  Bay. 
C  Comparison  q£   velocity  and  backscattering  data. 

Velocity  data  was  taken  at  the  two  locations  as 
described  in  Appendix  VI.  The  data  were  taken  with  the  intention 
of  comparing  with  velocity  data  at  the  same  time  utilizing  the 
computer  program  of  appendix  VII. 

Unfortunately  the  DuMont  404  B  pulse  generator  was  not 
available  when  taking  data  at  the  N.E.L.  tower  in  San  uiego.  A 
GR  1217-A  pulse  generator  was  used  instead.  Its  relatively 
unstable  pulse  width  caused  considerable  horizontal  jitter  reducing 
the  reading  accuracy  to  about  0.07  percent.  Figure  6=5?  shows 
velocity  data  and  simultaneous  data  derived  from  scattering  data 
taken  at  the  N.E.L.  tower.  The  broken  line  is  the  bubble  free 
velocity.  It  was  chosen  so  that  the  computed  dispersion  curve  best 
fits  the  data  points.  It  was  not  possible  to  determine  the  bubble 
free  velocity  because  simultaneous  salinity  readings  were  not  taken 
for  substitution  into  the  empirical  formula  for  sound  velocity  in 
sea  water,  [ll]  The  data  for  figure  6=5  were  taken  on  15  February 
196$  in  mid  afternoon  and  corresponds  to  figure  4-25  and  6  • 
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Attenuation  data  was  not  taken  simultaneously  with  velocity  data. 

Figure  6-6  is  a  dispersion  curve  computed  from 
scattering  data  taken  on  the  morning  of  18  March  1965  in  Monterey 
Bay.  The  points  appearing  on  Figure  6-6  are  velocity  data  points 
taken  at  the  same  time.  Again,  since  salinity  was  not  known  the 
bubble  free  velocity  can  not  be  determined.   Attenuation  data  as 
a  function  of  frequency  were  not  taken  in  Monterey  Bay. 

Examination  of  figures  6-*>  and  6-6  reveals  that  the 
velocity  data  do  not  disagree  with  the  velocity  dispersion  curves 
obtained  from  scattering  and  attenuation  data.  Unfortunately  the 
largest  velocity  shifts  in  these  dispersion  curves  occur  at 
frequencies  below  the  low  frequency  limit  of  the  velocity  measuring 
apparatus o  This  could  be  remedied  by  using  a  transducer  which  has 
a  flat  response  that  extends  to  lower  frequencies.   It  is  the 
transition  from  stiffness  to  resistance  control  (about  4.0  Kcps 
in  the  present  transducers)  that  produces  a  phase  shift  which 
prevents  accurate  measurement. 

Examination  of  figures  5-1,  6-f>  and  6-$'  reveals  that 
the  shape  of  the  dispersion  curves  obtained  from  the  Glotov  bubble 
distribution  and  the  various  bubble  distributions  obtained  in 
this  paper  are  very  ^similar.  Also  the  velocity  shift  from  bubble 
free  velocity  changes  very  little  with  changing  frequency  above 
80  Kcps.  If  we  assume  a  ^lotov  distribution  and  take  one  velocity 
measurement  at,  say,  80  Kcps  we  can  obtain  the  total  bubble  volume 
concentration  for  the  entire  distribution  by  substituting  the 
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measured  velocity,,  V  into  equation  6-1 

Vol.  Cone.  =9.71  (V-V0)  x  1CT8         (6-1) 
which  is  derived  from  figure  5-1  at  80  Kcps.  This  equation 
requires  a  knowledge  of  V0  the  bubble  free  velocity  which5 
according  to  figure  5-1^  should  be  obtained  below  three  Kcps0 
Unfortunately  the  instrumentation  described  in  this  paper  is 
not  operable  below  20  Kcps.  However s  this  single  measurement 
technique  could  be  useful  in  studying  the  decay  of  bubbles  in  a 
wake  with  time  by  simply  assuming  the  termianl  velocity 
measurement  is  the  bubble  free  velocity ^  V0. 
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7.  Conclusions. 

It  is  postulated  that  the  total  attenuation  measured  at  sea  at 
frequencies  below  60  Kcps  is  due  to  bubbles,  the  effect  of  the 
bubbles  increasing  as  frequency  is  reduced. 

At  all  frequencies  there  will  be  attenuation  with  negligible 
scattering  due  to  plankton  of  size  ka<0.1.  This  effect  is  dominant 
above  60  Kcps. 

It  is  impossible  to  assess  the  effect  of  any  suspended  inorganic 
material  in  the  sea  since  its  acoustic  effect  would  be  the  same  as 
that  of  plankton  (algae-type).  Since  the  San  ^iego  sea  water  was 
very  clear  it  is  assumed  that  there  was  a  negligible  amount  of 
suspended  particles. 

The  velocity  shift  measurements  do  not  disagree  with  the  above 
hypothesis. 
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APPENDIX  I 

The  Electrostatic  Transducer 
A.  Brief  Theory- 
Electrostatic  transducers  have  been  described  in  the  literature 
by  a  number  of  individuals. 

To  arrive  at  the  expected  forms  of  the  receiving  and  transmitting 
sensitivities  the  following  derivation  is  presented  starting  with 
the  equation  of  motion  for  any  point  localized  mechanical  system 
with  an  applied  periodic  force t 

md2x  ^  Rmdx  +  9X  =  p^jwt  (i_x) 

where  x  is  the  displacement. 
m  is  the  mass. 

Rjq  is  the  mechanical  resistance, 
s  is  the  stiffness 
F  is  the  force  applied. 
If  the  displacement  is  simple  harmonic  then: 

X  =  Xeswt 
and  (iV-l)  beomces: 

-Xw^  +  jXwRm  +  Xs  =  F  (1-2) 

Considering  the  transducer  as  a  receiver  a  substitution  can  be 
made  for  the  force: 
F  =  PA 
where  P  is  the  pressure 

A  is  the  area  of  transducer  face 
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and  further,  if  a  large  polarizing  voltage  is  applied  such 
that  the  total  charge  is  essentially  constant,  then  the  signal  voltage 
developed  (V)  is  approximately  proportional  to  the  displacement 
X: 

V  =  kX 
hence : 

-kVvAa  +  ,1  kVwRm  =  kVs  =  PA 

-Vw  fa   -  ra  (1-3) 

k 

Microphone  sensitivity  is  £iven  by: 

M  =  V 

P 

so; 

-Mw2m  +  jMwttm  '+  Ms  =  A  -  constant  (1-4) 

k 

We  have  thr<  t  indicate  three  distinct  and  different 

regions  of  operations 

(a)  Low  frequency  or  stiffness  controlled  regions 

Ms  5=«  constant 
M«;  Constant 

(b)  High  frequency  or  mass  controlled  region: 

Mw^m  -^  Constant 

M   Constant 
M^    W2 

In  this  region  the  response  should  fall  off  at  a  rate  of  12  db 

per  octave  with  increasing  frequency. 

(c)  Mid-frequency  or  resistance  controlled  regions 

MwRm  *=«  Constant 

M  .  Constant  (1-5) 
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In  this  region  the  hydrophone  response  should  fall  off  at  a 
rate  of  6  db  per  octave. 

The  E/S  transducer  with  its  large  diameter  to  wavelength  ratio 
is  assumed  to  be  operating  in  this  region  over  the  frequency  range 
of  the  investigation  (20  Kcps  to  200  Kcps). 

Considering  the  same  region,  one  way  of  expressing  the 
transmitting  sensitivity  is  J 

where  I  is  the  driving  current  for  constant  voltage  V, 
but  since: 

I  =  wCV 
where  C  is  the  capacitance  of  the  transducer 


S-  =  L. 

1   wCV 


but  from  (1-5)  Vpc — £-  (  constant  J 

W   v        ' 


so  Si  fe  Constant  =  constant 
G 

Another  way  of  expressing  the  transmitting  sensitivity  is; 

Sv  =  l 

where   V  is   the  driving  voltage  for  constant   current. 
Since 

V  =1. 
wC 

S     =  J^S         but  from   (1-5)    '/.      P (Constant) 
v  j 

w 


Sv     =  w  ^Constant)  (1-7)    £S  —■     f  constant) 
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In  this  case  the  transmitting  sensitivity  should  increase  with 
frequency  by  6  db  per  octave. 

The  calibration  data  obtained  for  5  transducers  is  in  agreement 
with  the  resistance-control  assumption  and  the  conclusions  of 
equations  (1-5),  (1-6)  and  (1-7). 

Looking  at  equations  (1-5  and  (1-6)  it  can  be  seen  that  the 

impedance  of  the  transducer  should  also  be  frequency  dependent: 

from  (IV-5)   M  =  •?   ^  constant 

w 


from  (IV-6)  S±  =  £  ^-Senstant 


Z  jti  ^OI)S  tant  (1-8) 

w 

This  is  to  be  expected  since  the  E/S  transducer  is  no  more 
than  a  capacitor. 
B,  Design  and  Construction 

1.  General 

Reference  is  made  to  Figure  1-1  which  shows  a  cross-sectional 
view  of  the  transducer.  The  fixed  electrode  is  the  aluminum  plate 
while  the  movable  electrode  is  the  thin  aluminum  film  on  the  outer 
side  of  the  mylar;  the  one  mil  mylar  constitutes  the  solid  dielectric. 
The  electrical  circuit  is  completed  by  means  of  a  coaxial  cable 
connected  to  the  two  electrodes  by  means  of  the  screw  terminal  and 
silver  painted  area  shown  on  the  diagram.  A  small  air  cushion  exists 
between  the  fixed  electrode  and  the  mylar;  the  concentric  grooves 
in  the  aluminum  plate  add  to  the  air  cushion  and  increase  the 
sensitivity  by  adding  compliance  to  the  system. 
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By  means  of  the  air  chamber  in  the  back  of  the  transducer 
pressure  equalization  is  obtained  which  makes  the  sensitivity 
independent  of  depth.  Such  pressure  equalization  was  described 
by  Kuttruff  and  Wille  (12  )j  there  will  of  course  be  a  limiting 
depth  for  any  given  size  of  chamber  but  one  of  the  transducers 
constructed  for  this  investigation  was  successfully  operated  at  a 
depth  of  40  meters  with  no  apparent  change  in  sensitivity. 

The  transducer  diameter  of  24  cm  is  equivalent  to  8  wavelengths 
at  50  kcps;  this  was  considered  satisfactory  and  made  for  a  con- 
venient handling  size. 

Many  hours  were  spend  in  developing  techniques  for  the 
construction  of  a  satisfactory  transducer;  so  many  subtle  details 
were  involved  in  the  final  technique  that  a  step  by  step  procedure 
is  outlined  in  Appendix  I-B.2. 

Once  the  mylar  was  bonded  to  the  plexiglass  rim  a  coating  was 
required  to  protect  the  thin  aluminum  film  from  the  action  of  the 
sea  water.  Such  a  coating  has  negligible  effect  on  the  transducer 
performance  provided  that  its  mass  is  not  significant,  its  thickness 
is  much  less  than  a  quarter  wavelength,  and  it  does  not  contain  air 
bubbles.  In  practive  a  coating  of  Epon  828  of  thickness  0.35  to 
Q.40  mm  was  consistently  obtained;  this  more  than  satisfied  the 
quarter  wavelength  criterion  even  at  the  upper  limit  frequency  of 
200  Kcps.  The  response  curves  obtained  for  several  transducers 
showed  that  the  added  mass  of  the  coating  did  not  effect  the 
resistance-control  assumption  described  earlier. 
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In  determining  the  resonant  frequency,  resistance  and 
stiffness  of  the  system  it  was  assumed: 
Rw  ^  S  at  20  Kcps 
R  =  Rm  +  Rr 
Rm  $   Rr 
R^te  Rr 
where  S  is  the  stiffness 

Rm  is  the  mechanical  resistance 
Rr  is  the  radiation  resistance 
for  ka»l  -*>   Rr  =  ^oCTTa2  =  2.05  x  10$ 


hence:  SfcsRrw  =  4.1  x  10^ 


Newton s 


Meter 
The  resonant  frequency  is  given  by: 


Wo 


/Mo  + 


Mr 

where  Mo  is  the  measured  mass  of  the   radiator  =  0.0228  kg. 
Mr  is  the  added  mass  due  to  the   radiation  reactance. 

let  Wo  =  -J =  4.24  x  106 

Mo 

For  a  piston  in  the  range  Ka^>  1,  Mr  ^   Mo 

hence  the  resonant  frequency  is : 
f o  a  680  Kcps 

Front  and  rear  views  of  a  finished  transducer  are  shown  in 
Figures  1-2  and  1-3  respectively. 

The  greatest  difficulty  encountered  was  in  maintaining  a 
water-tight  seal  for  the  face;  if  water  penetrated  the  protective 
epoxy  coating  it  would  attack  the  aluminum  film  and  gradually 
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reduce  the  effective;  area  of  the  transducer.  An  example  of  face 
deterioration  due  to  water  penetration  is  shown  in  Figure  I-45  a 
photograph  of  the  face  removed  from  transducer  4  after  about  a  week 
of  operation  in  the  sea. 

Two  different  resins  were  used  in  constructing  the  transducers, 
Shell  Epon  828  and  Shell  ^pon  1004.  Operational ly,  there  was  little 
difference  but  the  828   was  used  in  most  cases  because  of  its 
easier  handling  qualities.  Both  epons  are  descriDed  in  Appendix 
I-B.3. 

2.  Construction  Procedure 
Problem  Statements 

(1)  To  provide  sufficient  tension  in  the  mylar  face 
while  the  adhesive  is  drying, 

(2)  To  keep  water  from  getting  in  under  the  rubber 
pressure-compensating  diaphragm, 

(3)  To  seal  the  face  so  that  the  water  cannot  react 
with  the  aluminum  coating  on  the  mylar, 

(4)  To  keep  the  contained  air  from  being  forced  out 
through  the  coaxial  cable, 

(5)  To  provide  a  good  electrical  contact  between  the 
mylar  aluminum  film  and  the  outer  conductor  of  the  coaxial  cable • 

Construction  Procedures 

(1)  The  following  is  a  step  by  step  procedure  for  the 
final  stage  construction  of  the  transducer  using  Figure  1=1  diagram 
as  a  reference. 

(2)  The  area  around  the  stuffing  tube  hole  is  coated 
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with  silver  paint  and  allowed  to  dry. 

(3)  The  coaxial  cable  must  be  fitted  to  the  transducer 
by  means  of  the  stuffing  tube.  The  cable  is  fed  thru  the  stuffing 
tube,  the  shielding  layed  back  and  soldered  to  the  body  of  the 
stuffing  tube;  the  thread  of  the  stuffing  tube  is  coated  with  neoprene 
before  insertion?  the  portion  of  the  thread  that  protrudes  after 
tightening  must  be  left  bare  so  that  a  good  electrical  connection 

can  be  made  to  the  aluminum  film  on  the  face. 

(4)  After  the  stuffing  tube  connection  is  tightened, 
the  center  conductor  can  be  formed  into  a  loop  and  connected  to  the 
aluminum  plate  by  means  of  the  provided  screw. 

(5)  Since  the  "connection  cavity"  provides  a  path  for 
the  air  to  be  forced  out  it  is  filled  with  Epon;  this  effectively 
seals  off  the  contained  air. 

(6)  Epon  is  spread  on  the  back  of  the  transducer  on  the 
area  which  will  be  directly  under  the  diaphragm  claming  ring,;  this 
ensures  a  good  seal  when  the  diaphragm  is  in  place. 

(7)  A  sheet  of  rubber  is  placed  over  the  back  of  the 
transducer  and  the  clamping  ring  placed  in  position;  holes  are 
punched  thru  the  rubber  when  the  pre-drilled  holes  of  the  ring  are 
lined  up  with  the  holes  in  the  transducer  back;  the  rubber  is 
stretched  during  this  operation  and  held  taunt  while  the  ring  is 
tightened  down.  The  excess  rubber  can  be  trimmed  off  with  a  knife. 

(8)  The  ring  area  is  coated  with  Epon  to  complete  the 
sealing  of  the  back.  The  transducer  is  left  to  dry  21+   hours. 

(9)  A  piece  of  mylar  out  to  size  is  placed  on  the 
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aluminum  frame  with  the  aluminum  film  facing  up;  an  ohmeter  should 
be  used  to  determine  which  side  is  the  aluminum  side  since  it  is 
difficult  to  determine  this  visually. 

(10)  The  mylar  is  stretched  over  the  frame  by  means  of 
strips  of  electrical  plastic  tape  adhered  to  the  mylar  and  pulled 
down  and  over  the  edge  of  the  frame.   Care  should  be  taken  to  keep 
the  tape  off  that  part  of  the  mylar  within  the  area  of  the  frame  hole. 

(11 )  The  frame  is  put  aside  and  a  film  of  Epoxy  Glue 
applied  to  the  outer  ring  of  the  transducer  between  the  edge  and 
the  aluminum  plate  (or  edge  of  trough  if  later  model). 

(12)  The  frame  is  fitted  over  the  transducer  so  that  the 
mylar  rests  on  the  facej  the  edges  are  pressed  to  ensure  good  contact 
with  the  epoxy  glue^,  the  frame  is  then  weighted  to  ensure  tension 

in  the  mylar  while  the  glue  drys|  this  setup  should  be  left  overnight. 

(13)  The  mylar  can  be  roughly  cut  away  to  enable  removal 
of  the  framej  the  mylar  is  then  "cut  back"  from  the  edge  about  one 
eighth  of  an  inch  leaving  an  outer  ring  of  plexiglas. 

(14)  Silver  paint  is  not  applied  to  the  exposed  threads 
of  the  stuffing  tube  and  spread  down  over  the  front  edge  of  the 
transducer  and  over  a  small  area  of  the  aluminum  film0  Making  small 
groves  in  the  edge  at  the  area  where  the  paint  is  to  be  applied 
helps  in  providing  a  good  electrical  contact.  Once  the  silver  paint 
has  dried  a  continuity  check  should  be  made  with  an  ohmmeter. 

(15)  A  mixture  of  Epon  is  prepared  and  spread  over  the 

face  of  the  transducer s   which  has  been  placed  on  a  steady s   level 

surface,,  by  means  of  a  corrugated  cardboard  "spreader"  about  U'lx-5n 

(any  material  which  will  not  scratch  the  aluminum  film  will  suffice). 
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The  Epon  can  be  spread  over  the  entire  face  quite  uniformly  with 
minimum  thickness;  care  must  be  taken  to  ensure  that  there  are  no 
pinholes  exposing  the  aluminum  and  that  the  coating  is  spread  right 
to  the  edge  of  the  transducer. 

(16)  After  the  Epon  has  dried  the  capacitance  of  the 
transducer  can  be  measured  on  a  bridge;  it  should  be  ^  0.005  uf. 

(17)  The  transducer  is  now  ready  for  use. 
3.  Epon  Characteristics 

Epons  828  and  1004  were  both  used  in  transducer  construction 
with  curing  Agent  U;  both  were  satisfactory  but  the  Epon  828 
being  a  liquid  was  easier  to  use  and  became  the  perferred  coating. 
The  characteristics  of  the  Epons  are  as  follows  1 

Epon  828:  A  pourable  liquid  at  room  temperature;  widely 
used  in  wet  lay-up  laminating^  casting  and  the 
pathing  of  electrical  equipment ,  and  in  form- 
ulating adhesives  for  bonding  wooda  plastics,, 
metal  and  glass  to  themselves  and  to  each  others 
its  low  viscosity  makes  it  attractive  for 
contact  pressure  and  vacuum  bag  laminating  at 
room  temperature. 
Epon  1104s  Supplied  as  a  solid  resin  and  requires  a  solvent 
such  as  Acetone.  Combines  flexibility^,  adhesion 
and  toughness  with  good  chemical  resistance 
properties. 
Curing  Agent  Us  Provides  fast  cure  at  room  temperature;  has 
short  pot  life  and  low  toxicity;  gives  strong 

lamination. 
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Mixtures   (by  weight):  828 _  5, 

Agent  U         1 

lQOii  1£ 

Agent  U         1 

C.   Characteristics 

1.  Response 

With  the  help  of  the  NEL  Transducer  Evaluation  Center 
(TRANSDEC)  calibration  data  was  obtained  for  5  of  the  6  transducers 
taken  to  San  Diego  for  field  measurements. 

The  equations  developed  in  Appendix  I-A  indicated  that  the 
transmitting  sensitivity  with  respect  to  driving  current s   in 
microbars  of  pressure  at  one  meter  per  ampers  input  should  be  a 
constant  in  the  resistance  controlled  regionj  Figure  1-5  shows 
such  curves  for  the  5  transducers  calibrated  and  with  the  exception 
of  transducer  6  are  within  1.5  db  of  the  89  db  level  throughout 
the  20  Kcps  to  150  Kcps  frequency  range  of  the  calibration. 

The  transmitting  sensitivity  can  also  be  expressed  in  terms  of 
microbars  per  volt  and  as  was  shown  in  Appendix  I-A  the  sensitivity 
should  vary  directly  as  the  frequency  or  increase  at  a  rate  of  6 
db  per  octave |  Figure  1-5  shows  that  such  is  the  behavior  for  all 
6  transducers. 

The  receiving  response  should,,  from  the  developed  theory <,  vary 
inversely  as  the  frequency  or  fall  at  a  rate  of  6  db  per  octave. 
Figure  1-7  shows  that  the  6  db  per  octave  behavior  is  followed 
reasonably  well  although  the  fall  rates  of  transducers  5X  and  6  is 
slightly  different  from  that  of  the  other  3  transducers  in  that 
the  slopes  of  the  response  curves  are  slightly  different.  This 
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must  be  attributable  to  the  fact  that  transducers  5^  and  6  were  faced 
with  Epon  828  while  1,  2  and  3  were  faced  with  Epon  1004. 

2.  Radiation  Characteristics 

The  fact  that  the  transducer  diameter  is  quite  large  (  8 
wavelengths  at  50  Kcps)  makes  it  highly  directive.  The  directivity 
patterns  for  the  5  transducers  calibrated  are  shown  in  Figures  1-12 
to  1-18;  for  all  5  were  calibrated  at  60  Kcps  and  in  addition 
patterns  for  transducer  3  were  obtained  at  100  Kcps  and  150  Kcps. 

In  comparing  the  radiation  behavior  of  the  transducers  with 
that  of  the  ideal  flat  piston  three  things  were  considered:  the 
effective  radius,  the  db  differences  between  the  main  lobe  and  the 
first  side  lobes,  and  the  directivity  index. 

The  calculated  effective  piston  radii  ranged  from  12.5  to 
15.5  cm  with  a  mean  of  13.2  cmj  this  was  slightly  greater  than  the 
actual  radius  of  12  cm  and  is  not  surprising  since  the  face  must  be 
subject  to  some  flexing  action  which  results  in  a  pattern  normally 
expected  from  a  larger  piston. 

The  directivity  index  ranged  from  29.7  db  to  31. 8  db  with  a 
mean  of  30.5  dbj  the  theoretical  value  for  a  piston  is  29.6  db. 

The  db  difference  between  the  main  lobe  and  the  first  side 
lobes  proved  to  be  much  more  inconsistent.  The  differences  ranged 
from  10.6  db  to  16.6  db  with  a  mean  of  13.8  db.  The  theoretical 
difference  for  a  piston  is  17.6  db.  In  seeking  an  explanation  it 
would  seem  that  the  fact  that  two  different  epons  were  used  in  the 
face  construction  would  be  a  factor,  but  the  two  transducers 
closest  to  the  theoretical  value  (16.6  db  and  16.2  db)  had  different 
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faces,  so  some  other  non-obvious  feature  of  the  construction  must 
be  at  fiiult;  the  thickness  of  the  epon  at  the  edges  which  varied 
from  the  mean  thickness  by  as  much  at  15$  might  be  a  significant 
factor.   The  variation  in  the  epon  thickness  and  the  fact  that  the 
epoxy  glue  was  sometimes  squeezed  in  over  the  edge  of  the  aluminum 
plate  thereby  distorting  the  circular  shape  of  the  radiator  slightly 
accounts  for  most  of  the  assymmetry  in  the  radiation  patterns. 

In  use  all  the  back  radiation  portions  of  the  directivity 
patterns  were  probably  suppressed  by  the  rigid  steel  mounting  plate. 
The  mounting  can  be  seen  in  Figure  1-8. 

3.  Capacitance  and  Impedance  Characteristics. 

The  electrostatic  transducer,  like  the  condenser  microphone, 
is  a  device  that  depends  for  its  operation  on  the  value  of 
capacitance  between  a  fixed  plate  and  a  stretched  membrane  which 
in  this  case  is  the  aluminum  film  on  the  mylar. 

As  in  any  capacitors 

C  oC  l/d  where  d  is  the  separation  between  plates.  For 
the  E/S  transducer  the  dielectric  is  made  up  of  the  mylar  and  the 
air  custion  in  series  so  that  the  capacitance  C  is  given  by: 

1  =  1-  +  jL 

where  C^  is  the  capacitance  due  to  the  mylar  part  of  the 
dielectric. 
C^  is  the  capacitance  due  to  the  air  part  of  the 
dielectric. 
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finally 

C  =  J%&A (1-9) 

'vd^Km  +  dm) 

where  Km  is  the  dielectric  constant  for  mylar  =  3 

fio  is  the  permittivity  of  free  space  =  8.85  x  lcH-2 
A  is  the  area  of  the  aluminum  plate  =  4.53  x  1CP "  m^ 
d^  is  the  thickness  of  the  air  custion 
dy  is  the  thickness  of  the  mylar  =  2.54  x  lO™-8  m  (l  mil) 

Using  a  measured  value  of  capacitance  at  a  depth  of  10  feet 
it  was  found  that  the  thickness  of  the  air  cushion  at  that  depth 
was  1.88  x  1Ct5  m  (0.741  mil). 

It  was  found  that  the  value  of  capacitance  increased  with 
water  depth  due  to  compression  of  the  air  cushiony  the  increase 
was  quite  rapid  to  5  feet  and  then  leveled  off  to  a  linear  increase 
of  approximately  100  uuf  per  foot  of  depth.  If  the  pressure 
compensation  system  were  perfect  the  capacitance  would  not  change 
at  all  with  depth  but  such  a  system  would  require  an  exact  combination 
of  back  diaphragm  and  mylar  membrane  stiffness  relationships. 

Capacitance  data  taken  at  the  NEL  Tower  to  a  depth  of  45  feet 
is  plotted  in  Figure  1=9.  Curve  data  between  the  surface  and  5 
feet  is  only  approximate  due  to  difficulties  created  by  wave  action. 

Shallow-water  data  was  later  gathered  in  the  USNPGS  anechoic 
tank;  this  data  is  plotted  in  Figure  1-10. 
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The  impedance  of  the  E/S  transducer  depends  solely  on  the 
value  of  capacitance  discussed  above.  A  curve  showing  impedance 
variation  over  the  frequency  range  of  interest  is  plotted  in 
Figure  1-11  using  the  average  value  of  capacitance  at  10  feet  from 
Figure  1-9. 
D.  Factors  Affecting  Sensitivity 

The  receiving  sensitivity  of  the  E/S  transducer  depends  on  a 
number  of  factors  which  can  be  expressed  as  follows: 

E„a  E 
Sensitivity  »M«   C-  (1-10) 

where  Eq  is  the  dc  polarizing  voltage 
ae  is  the  effective  radius 
E  is  the  source  pulse  amplitude 
d  is  the  spacing  between  the  aluminum  plate  and  the 

aluminum  film. 
T  is  the  tension  in  the  stretched  mylar. 
Zjf  is  the  mechanical  input  impedance, 
e  is  the  charging  factor  associated  with  the  mylar 
The  polarizing  voltage  Eq  was  maintained  at  450  volts  (upper 
limit  of  available  power  supply)  throughout  the  investigation  while 
the  source  pulse  voltage  (E)  used  was  50  volts  peak;  these  values 
gave  an  approximate  10 si  ratio  and  ensured  that  harmonic  distortion 
would  be  minimized. 

The  ratio  Eg/d,  the  electric  field  intensity,  will  increase 
with  depth  as  the  air  cushion  is  compressed  as  described  earlier; 
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such  action  increases  the  sensitivity. 

In  construction  of  the  transducer  pre-stretching  the  mylar 
determined  the  value  of  the  tension  Tj  this  was  a  compromise 
situation  since  the  stiffer  the  mylar  was  the  more  the  transducer 
behaved  like  an  ideal  flat  piston  and  it  is  also  a  factor  in  the 
pressure  compensation  system;  on  the  other  hand  increasing  T 
decreases  the  sensitivity  as  shown  in  (1-10), 

Increasing  the  depth  increases  the  mechanical  input  impedance 
Z^  which  tends  to  decrease  the  sensitivity,  however  the  increase 
of  electric  field  intensity  with  depth  is  the  dominant  factor  so 
that  the  overall  effect  of  increased  depth  is  an  increase  in 
sensitivity,  disregarding  the  pressure-equalization  system. 

It  was  found  that  there  was  a  charging  factor  (e)  associated 
with  the  mylar  due  to  the  polarizing  voltage  Eq.  For  a  given 
polarity  of  Eq  for  a  period  of  time  would  have  a  neutralizing 
effect  so  that  if  the  polarity  was  again  reversed  there  was  an 
observable  increase  in  sensitivity.  This  effect  did  not  present 
any  major  problem  and  did  not  have  any  noticeable  effect  on  the  data 
as  long  as  the  polarity  remained  unchanged  for  the  duration  of  a 
data  taking  period. 

In  addition  to  all  of  the  above,  the  fact  that  the  capacitance 
is  depth  sensitive  as  described  earlier,  means  that  wave  action 
could  also  be  a  factor  affecting  sensitivity,  particularly  if  the 
transducer  is  operated  close  to  the  surface.  As  an  example  of  this 
effect  consider  the  following? 
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Wave  velocity  =  10  knots 
Wave  period  =  3  sec 

Wave  amplitude  =  6  f t  (crest  to  trough) 

For  a  wave  train  equivalent  to  20  msec  the  maximum  change  in 
depth  for  the  period  of  the  wavetrain  works  out  to  be  ~  0.1  ft; 
the  resulting  sensitivity  changes  are  tabulated  in  Table  1-1. 

Table  1-1 
Depth-ft       C-uuf         C-uuf         M 


5 

15400 

45 

.0260  db 

10 

16400 

10 

.0054  db 

20 

17500 

10 

.0050  db 

40 

19600 

10 

.0044  db 

Since  the  transducer  was  always  operated  at  depths  of  10  feet 
or  greater  the  above  effect  was  not  significant  in  this  investigation, 
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1  db  per  radial  division 
Effective   Radius     12.5  cm  Frequency  60  KCS 

1st  Side   Lobes  Down     13.5  db  Test  Distance   8  meters 

Directivity  Index     29.9  db  Depth  6.1  meters 

Figure    1-12  Directivity  Pattern-Transducer  No.    1 
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Effective  Radius   15.5  cm 
1st  Side  Lobes  Down  10,6  db 
Directivity  Index  31.8  db 


Frequency  60  KC3 

Test  Distance  8  meters 

Depth  6  meters 


Figure   I-l3     Directivity  Pattern-Transducer  No,   2 
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1  db  per  radial  division 


Effective   Radius   12.3  cm 
1st   Side   Lobes  Down  16.6  db 
Directivity  Index  29.7  db 


Frequency  60  KCS 

Test  Distance   6  meters 

Depth  6.1  meters 


Figure   1-14     Directivity  Pattern-Transducer  No.   3 
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1  db  per  radial  division 

■ 

Frequency  100  KCS 
Test  Distance   8  meters 
Depth  6.1  meters 

Figure   1-15     Directivity  Pattern-Transducer  No.   3 


143 


t  db  per  radial  division 

Frequency  150  KCS 
Test  Distance   8  meters 
Depth  6.1  meters 

Figure   1-16     Directivity  Pattern-Transducer  No*   3 
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1  db  per  radial  division 


Effective  Radius  12.3  cm 
1st  Side  Lobes  Down  16.2  db 
Directivity  Index  29.7  db 


Frequency  60  KCS 

Test  Distance  8  meters 

Depth  6.1  meters 


Figure  1-17  Directivity  Pat tern- Transducer  No.  5X 
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1  db  per  radial  division 


Effective  Radius  13.5  cm 
1st  Side  Lobes  Down  12,2  db 
Directivity  Index  30.6  db 


Frequency  60  KCS 

Test  Distance  8  meters 

Depth  6.1  meters 


Figure  1-18  Directivity  Pat tern- Transducer  No.. 6 
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APPENDIX  II 
SIGNAL  GATE 

This  circuit  utilizes  a  General  Radio  1217-B  Pulse  Generator 
to  gate  an  AC  input  signal.  When  used  in  conjunction  with  the 
General  Radio  13V6-A  Tone  Burst  Generator^,  it  provides  an 
additional  50  db  separation  of  on  and  off  signals  and  it  provides 
more  flexibility  in  controlling  pulse  width.  The  circuit  is 
shown  in  Figure  II -1. 

The  gate  consists  of  two  sections.  The  first  section  is  a 

common  emitter  amplifier  whose  collector  bias  voltage  is  supplied 

by  the  negative  terminals  of  the  pulse  generator.  The  width 

of  the  collector  bias  pulse  determines  the  width  of  the  resulting 

wave  pulse  output.  The  bias  resistors  are  adjusted  for  maximum 

gain.  Since  high  gain  is  the  principal  means  for  obtaining  a  high 

on/off  signal  ratio  a  high  gain  transistor  was  used  for  this 

amplifier.  The  forward  DG  current  transfer  ratio  for  the  2N1380 

is  listed  as 

30  minimum 
hFE  =  100  typical 
300  maximum 

The  purpose  of  G^_  is  to  present  a  high  impedance  to  the  signal 

generator  and  provide  DG  isolation.  The  value  of  Cj_  was  determined 

by  decreasing  it  until  the  peaks  of  the  first  cycle  of  a  200  Kcps 

wave  pulse  were  the  same  amplitude  as  the  subsequent  peaks. 

The  collector  output  voltage  contains  a  negative  five  volt 

pedestal  which  is  removed  by  the  second  section 0  G3  and  D-j_  clamp 
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the  "plus"  output  of  the  General  Radio  1217-B  pulse  generator  so 
that  the  minimum  value  of  the  pulse  is  zero  at  the  cathode  of  Dj_. 
This  positive  going  pulse  is  added  to  the  wave  pulse  thereby 
nullifying  the  negative  pedestal.  R5  provides  a  fine  control  of  the 
pedestal  height.  The  addition  is  accomplished  by  the  resistive 
network  consisting  of  R^,  R5,  R^,  and  Ry.  Figure  II-2  shows 
typical  wave  pulse  outputs  at  40  and  200  Kcps. 

For  optimum  performance  the  collector  bias  pulse  amplitude 
should  be  adjusted  so  that  the  gated  signal  is  clipped  symmetrically 
when  the  amplifier  is  overdriven.  Then  the  oscillator  amplitude 
should  be  adjusted  for  maximum  output  amplitude  without  clipping. 
These  adjustments  should  be  made  every  time  the  frequency  is  changed. 
When  time  does  not  permit  following  this  procedure,  the  pulse 
generator  should  be  set  once  at  ten  volts  and  the  oscillator 
amplitude  set  as  shown  in  Figure  II-3.  This  shorter  procedure  will 
cost  a  few  db  in  on/off  signal  ratio. 

Brief  specifications  are  as  follows; 

On/off  ration  =  61  db  at  25  Kcps 
56  db  at  60  Kcps 
45  db  at  200  Kcps 

Minimum  pulse  width  =  5  microseconds 

Maximum  pulse  width  =  8  milliseconds 
Figure  II-4  shows  the  variation  of  voltage  transfer  ratio  with 
frequency. 
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APPENDIX  III 
TRANSMIT-RECEIVE  SWITCHING  SECTION 

This  section  isolates  the  transmitter  from  the  reciver  using 
three  networks.  See  figure  III-l. 

The  network  consisting  of  D^  and  D£  offers  a  low  series 
impedance  to  high  voltage  signals  and  a  high  impedance  to  low 
voltage  signals.  Therefore  the  echo  signals  do  not  see  the  low 
power  amplifier  impedance.  It?  and  C]_  isolate  the  polarizing 
voltage  supply  from  the  AC  signal  voltages.  The  network  consisting 
of  Do  and  D,  presents  a  low  impedance  to  ground  for  high  voltage 
signals  and  a  high  impedance  to  ground  for  low  voltage  signals.  The 
main  bang  therefore  appears  across  R^0  Several  volts  of  main  bang 
to  appear  across  D3  and  D^  which  seriously  saturates  the  filter, 
amplifiers,  and  oscilloscope  display  thus  making  the  electronic 
switch  necessary.  The  network  at  the  output  of  the  GR-1217-A  pulse 
generator  differentiates  both  pulse  generator  outputs  and  clips 
the  resulting  positive  pulses.  The  remaining  negative  pulses 
trigger  the  electronic  switch  so  that  channel  A  is  on  during  the 
main  bang  and  channel  B  is  on  at  all  other  times.  The  switching 
transients  appearing  at  the  output  of  the  electronic  switch  are 
shown  in  Figure  2-4.  They  were  not  wide  enough  to  saturate  the 
receiver.  The  transmit-receive  switching  section  provides  about 
60  db  of  rejection. 
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APPENDIX  IV 
SYSTEM  CHAKACTEHISTIGS 

Overall  characteristics  of  the  pulse  echo  system  as  one 
unit  are  presented  below. 

The  usable  frequency  range  is  16  Kcps  to  200  Kcps.  The  upper 
frequency  is  limited  by  the  band  pass  filter,  and  low  electrical 
signal  to  noise  ratio..  Power  amplifier  distortion  at  high 
frequencies  was  noticable  but  not  serious.  The  low  frequency  is 
limited  primarily  by  the  divergence  of  the  transducer  beam. 

The  received  signal  to  noise  ratio  as  a  function  of  frequency 
is  shown  in  Figure  17-1.  Although  the  broken  curve  is  a  true 
indication  of  receiver  performance,  the  solid  curve  is  more  an 
indication  of  system  usefulness  since  the  fourth  echo  was 
arbitrarily  fixed  as  the  first  usable  echo  and  it  includes  transducer 
divergence  effects.  The  electronic  switch  was  the  primary  offender 
in  contributing  noise  to  the  system.  The  transmitted  signal  to 
noise  ratio  is  greater  than  90  db  below  60  Kcps  and  greater  than 
80  db  above  60  Kcps. 

The  maximum  wave  pulse  width  available  for  the  system  is  limited 
by  the  tone  burst  generator  to  128  cycles.  The  pulse  repetition  rate 
is  controlled  from  the  tone  burst  generator  and  may  be  varied  from 
a  continuous  signal  to  one  pulse  in  ten  seconds.  The  minimum  pulse 
width  is  five  microseconds  and  is  limited  by  the  signal  gate. 

The  receiver  output  deviates  from  linearity  no  more  than  about 
0.2$  db  or  2.8  percent.  This  figure  is  approximately  the  order  of  the 
resolution  of  the  instruments  used  in  making  the  linearity  check, 
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APPENDIX  V 
EQUIPMENT  SPECIFICATIONS 

The  following  is  a  list  of  specifications  of  the  commercial 
equipment  used  in  this  project.  Only  those  features  pertaining 
to  the  project  are  included. 
Audio  Oscillator  TS-382-DU 

Manufactured  by  TRAV-LER  Radio  Corporation 

Serial  number:  1152 

Frequency  ranges  20-200,000  CPS,  continuously  variable 
within  four  decade  ranges 

Output  signal  amplitude;  0-10  volts,  continuously  variable 
within  seven  decade  ranges 

Harmonic  distortion:  less  than  one  percent 

Output  impedance :  1000  ohms 
Tone  Burst  Generator 

general  Radio  Company  1396-A 

Serial  Number:  397 

Signal  frequency  range:  DC-500  Kcps 

Maximum  input  voltages  +  seven  volts 

Input  impedance:  About  10  kilohms 

Maximum  voltage  output:  +  seven  volts 

Total  distortion:  less  than  -60  db  at  1  Kcps  and  10  Kcps 

On/off  ratio:  -40  db  with  maximum  signal  input 

Output  impedance:  600  ohms 

Gating  voltage  output:  12  volts,  gate  open 

-12  volts,  gate  closed 
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Pulse  Generator 

General  Radio  Company  1217-B 

Pulse  repetition  frequency;   2.5  GPS  to  500  Kcps 

internally  generated.  Externally  controlled^  DC  to  1 
MCPS.  Output  pulse  is  started  by  negative  slope „ 

Output  pulse  durations  100  nsec  to  1  second  in  seven  decade 
ranges 

Rise  and  fall  times;  Less  than  20  nsec. 

Overshoot  and  noise;  less  than  5  percent  of  amplitude. 

PRF  stability;   .05  percent 

Pulse  duration  stability;   .03  percent 
Power  Amplifier 

General  Radio  Company  1233-A 

Serial  number;  539 

Frequency  Ranges   20  to  3  MC 

Output  power;  15  watts 

Voltage  output;   .1  microsecond  rise  time^,  negligible  overshoot. 

Input  voltage;   .2  volt  for  full  output 

Input  impedances  100  K-ohms5  37  pf 

Distortion;  Less  than  three  percent  at  rated  output. 

Noises  70  db  below  15  watts  or  .3  volt  peak  to  peak. 

Power  consumption;  120  watts  at  zero  output 
Electronic  Switch 

DuMont  Laboratories  type  330 

Serial  number;  1111 

Input  impedance;  2  megohms s   45  pf 
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Response:  Flat  to  DC  with  60  pf  load,  response  down  no  more 

than  1  db  at  10  MC. 
Hum  and  noise  levels  2.$  MV  maximum  exclusive  of  s  switching 

transients. 
Switching  signal  rise  time:  Two  percent  of  square  wave  period 
External  trigger  switching  rate:  Zero  to  100  Kcps 
Synchronizing  Voltage:  At  least  +  2  volts  peark  with  rise. 

.02  V/  s. 
Power  consumed:  160  watts 
Electronic  Filter 

Spencer-Kennedy  model  302 

Serial  Number:  1483 

Frequency  limits:  20  GPS  to  200  Kcps  in  four  decade  ranges 

Kate  of  attenuation  in  rejection  band:  18  db/octave  per  section 

Maximum  attenuation:  90  db 

Input  impedance:  2  megohms,  30  pf 

Maximum  input  signal:  Six  volts  peak 

Output  impedance:  300  ohms 

Minimum  external  load  for  operation  according  to  specification: 

Ten  K-ohms 
Noise  Level:  80  db  below  one  volt  (using,  detector  bandwidth 

from  zero  to  40  Kcps) 
Power  consumed:  120  watts 
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Amplifier 

Hewlett-Packard  model  450-A 

Serial  numbers:   228-06357  and  228-06379 

Gain:  20  db  +  1/8  db  or  40  db  +  1/8  db 

Frequency  responses  20  db  gains  within  -  g  db  5  to  1  MC 

40  db  gain:  within±5  db  to  10  to  1  MC 
Maximum  input  voltage:  20  db  gain:  1  volt 

40  db  gains  .1  volt 
Input  impedances   one  megohm^  15  pf . 
Outputs  Tenvolts  maximum  across  3000  ohms  (or  more) 
Distortion:  Less  than  one  percent  at  rated  output 
Internal  impedances  less  than  150  ohms 
Equivalent  input  noise  levels  40  db  gains  about  40   v 

20  db  gains  about  250  v 

Power  consumed s  55  watts 
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APPENDIX  VI 


TEST  SITES 


a)  The  U.  S.  Navy  Electronics  Laboratory  Oceanographic  Research  Tower. 

The  Tower  was  chosen  as  the  site  for  field  work  because  it 
provided  a  stable  platform  J29j  for  the  experimental  equipment  and 
was  situated  in  water  of  known  characteristics. 

The  Tower  is  located  one  mile  off  Mission  Beach  where  the 
bottom  is  gently  sloping  rippled  sand  and  the  mean  tidal  depth  is 
56  feet.  The  dominant  longshore  current  is  from  the  north  [_3Q] 
so  that  water  from  Mission  Bay  or  the  flood  control  channel  is 
unlikely  to  reach  the  Tower  and  the  surrounding  water  is  of  almost 
constant  salinity  in  the  winter. 

The  electronic  equipment  was  placed  on  the  Main  Deck  in  the 
room  at  the  right  hand  side  of  Figure  VI-1  with  the  transducers, 
reflector,  temperature  and  oxygen  sensors  mounted  on  a  car  beneath 
Figure  VI-2.  This  car  was  firmly  mounted  on  railway  tracks  and  could 
be  lowered  to  the  sea  bottom  so  that  it  provided  a  stable  platform 
at  any  desired  depth. 

The  low  self -noise  of  the  Tower  was  an  advantage  which  may  have 
been  somewhate  offset  by  echoes  from  the  structural  members  and 
marine  growth  shown  in  Figure  VI-2. 

Full  details  of  the  construction  and  operation  of  the  Tower 
will  be  given  in  a  publication  which  is  now  being  prepared  by  Dr. 
E.  C.  LaFond  of  NEL. 
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b)  Monterey  Bay 

Experiments  were  performed  in  Monterey  Bay  using  the  USNPGS 
Oceano graphic  Research  Vessel  as  platform.  Engines  were  stopped 
while  measurements  were  being  made  so  that  effects  due  to  wake 
and  screw  action  were  minimised.  The  vessel  was  usually  moored 
300  yards  off  shore  and  300  yards  east  of  Monterey  Municipal 
Wharf  in  27  feet  of  water  over  a  sandy  bottom. 
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APPENDIX  VII 
VELOCITY  SHIFT  COMPUTER  PROGRAM 

This  computer  program  computes  sound  velocity  shift  as  a 
function  of  frequency  in  sea  water  containing  a  distribution  of 
bubble  sizes.  The  basis  of  the  computation  is  the  expression  for 
added  compressibility  due  to  bubbles 


^  'efc^TO.*.  7K  r.co.l^       (5-4) 


ISC-']  *[<*#] 


For  a  distribution  of  discrete  bubble  sizes  the  compressibilities 
will  add   obtaining  at  frequency  ft) . 


•■'M-2.22..2L. 


i  rtf*.  Lsz  -ijSitgij* 


where  i  will  be  the  index  for  each  bubble  size.  The  velocity  at 
frequency,  W,  will  then  be; 


v(to)  = 


-^  +  a  (u) 

Co*" 


The  damping  factor  as  a  function  of  bubble  resonant  frequency 
has  been  approximated  by 


fx~103  J 
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di  "  W        gTT-"'  I  (VII-D 


Equation  VII-1  fits  Devin's  curve  within  about  five  percent  below 
200  Kcps.  03  . 

The  bubble  resonant  frequency  as  a  function  of  bubble  radius 
at  one  atmosphere  has  been  approximated  by 

8TT  x  106  (m_2) 

1     H 

where  Rj_  is  in  microns. 

In  equation  VII-2  the  bubble  is  assumed  to  behave  adiabatically, 
surface  tension  is  ignored  and  the  external  pressure  is  taken  as 
standard  sea  level  pressure. 

The  computer  program  provides  a  print  out  of  volume  concentration, 
formula  coefficient,  damping  factor,  and  resonant  frequency  for 
each  bubble  sizej  and  a  print  out  of  sound  velocity  shift  from 
1500  M/S,  density  times  compressibility  and  frequency.  The  formula 
coefficient  is 

c4  -  3  \/X±**02 

The  U.  S.  Naval  Postgraduate  School  subroutine  DRAW  was  utilized 
providing  a  graph  of  velocity  shift  versus  frequency.  [31J 
Figure  VII-1  is  a  flow  diagram  for  this  program.  Table  VII-1  is 
a  typical  print  out. 

Several  volume  distributions  were  put  into  the  program.  Each 
will  be  described  below. 
Glotov,  Kolobaev  and  Neuimin  distribution. 
For  25<-R<69  m. 

Vi  =  B(.2067Ri  -  5.17)  x  4.19  x  10-18R? 
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For  71^R^105/^ 

?i  =  B  x  1.91  x  lO"13^ 
For  107<R<325/m. 

V±   =  B  x  2.57  x  10~7A- 
where 

B  =  scaling  factor 

Vj  =  volume  concentration  of  all  bubbles  with  radii  one  micron 
on  either  side  of  Rj_. 
The  above  formulas  are  derived  from  information  given  by  Glotov  et. 

al.  in  their  paper.  However  the  authors  indicate  that  R^  should 

p 

appear  in  the  denominator  of  the  third  formula.  R^  was  used 

instead  because  Rj_  does  not  cause  V^  to  converge  rapidly  enough 
causing  a  large  discontinuity  at  the  larget  radius  of  325  microns. 
Distributions  from  attenuation  data. 

The  folume  concentrations  derived  from  figure  3-12  were 
normalized  to  two  micron  increments.  Two  distributions  were 
extracted  from  this  data;  the  first  utilizes  all  the  data  and  the 
second  utilizes  the  center  portion  of  the  distribution  only. 
First  curve: 
For  16<R<26  /^ 

V±  =  (,50Ri  -  8.0)  x  HT10 
for  28<R<60^ 

V±  =  U.912%  /  30.1)  x  10-10 
for  62«-R<-84/v 

V±  =  (,259Ri  -  15.05)  x  10-10 
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for     86<R«:110/v 

V±  =  4.3  x  10-VRi3 

for  112*.  R<  156  ^ 

V±  =  (O.IRi  -  7.8)  x  10-10 

for  158<-R<416/^ 

V±  =  3.05  x  10-3/1^3 

Second  Curve: 

For     58<.R<84//^ 


V±  =  (.259Ri  -  15.05)  x 


10' 


-10 


for     86-<.R<.458/a^ 

V±  =  4.3  x  HrtyR^ 
Distribution  from  Scattering  Data 

The  following  formulas  apply  to  Figure  6-2  and  6-4. 
San  Diego 
For  65<-R<149/^ 

V±  =  (.283Ri  -  18.3)  x  10"10 
for    151<*  R-*  465/^ 
V±  =  0.0080/Ri3 
Monterey  Bay 
For     65<R*-135/^ 

7±  -  (O.IRi  -  6.5)  x  10-10 
for  137*-R<-465/v 

Vi  =  0.0017/Ri3 
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42)21,21,2 
)=(.283*R( 

30 
)=.0080/R( 
3.0*VGL( I ) 

SOI  ,  (  I,  H 
T(  I5,lPi4El 
)=R(  D+2.0 
NUE 

0/22500C0. 
6283.2 

J=1,200 
J)=0.0 
0.0 

1  1=1, 2C0 
C(  I)«(  (WR( 
)«WR(I )/W( 
SUMV+VCL( I 
J)=SUMA(J) 
(  1.0/SGRTF 
=  M  J)/6283 
)=W( J)+628 
NUE 

400,  (SUV 
T130H0TCTA 

500 
T(3CH1VELC 

501 
T(26H     I 

600 
T(54H  FREG 

700, (EE( J 
T(0PF8.2,6 

10CC 
T(  1H1///) 

600 

700,(EE( J 

1000 

600 

700, (EEU 


IMPUTES    VELOCITY    SHIFT 
NEIL»S    SCATTERING    DATA 

1),R(201)!C(20  1),C(2  01),W(2C1),A(201), 
),  ITITLE(l2),EE(2Cl)  ,  VOL  (201  ) 


RESCN/NT  FREQ.   DAMPING  FACTOR   COEFFICIENT 


ACIUS 

RESONANT  FREGUENCY  (RACIANS) 

AMPING  FACTOR 

COEFFICIENT 

CONCENTPATICN  FCR  BUBBLE  SIZE  I 

INCEX 
Y  AS  INDEPENCANT  VARIABLE  (RACIANS) 
X  ACDEC  COMPRESSIBILITY  CUE  TC  BUBBLES  OF  SIZE  I 
CCMPRESS  1BILITY  CF  WATER  WITHOUT  BUBBLES 
F  ACDEC  COMPRESSIBILITIES  OVER  ALL  BUBBLES  SIZES 
LENCY  W(J) 
LOCITY  SHIFT  FROM.  1500  M/S  AT  FREGUENCY  W(J) 

VARIABLE  INCEX 
CY  IN  KCFS 


6«1000COC.C/R( I ) 
I)/6283.2)««0.322 
2 
I)-18.3)«1.0E-1G 


I)«»3 

/( (R( I)*1.GE-6«WR( I) )*»2) 

R( I  ), D( I  ),C(  I)»VCL(I ) ) 

6.7) 


I)/W(  J))«»2-1.)/(((WR(I)/W(J))««2-1.)»»2 

J)  )»«2) 

) 

+  A(  I) 

(AA  +  SUMAU)) J-15CC.0 

.2 

1.2 

V) 

L  VOLUME  CONCENTRATION  =  1PE16.7) 

CITY  SHIFT  DUE  TC  BUBBLES/) 

N  MCNEIL  DISTRIBUTION//) 

UENCY       RHO  X  COMPRESSIBILITY 


VELOCITY  SHIFT 


),SUMA(J),V( J),J=1,67) 

H   KC    1PE17.7.CPF18.7,  5H   M/S) 


),SUMA(J),V< J),J=68, 134) 

),SUMA(J),V< J),J=135,200) 
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ITITLEM  )=PHSCUNC 

ITITLE(2)=8HVELCCITY 

ITITLE(3)=8h  IN 

ITITLE(U)=8HWATER 

ITITLE(5)=0H  WITH 

ITITLE(6)=8HBUeeLES 

ITITLE(7)=8H 

ITITLE(8)=eHVERSUS 

ITITLE(9)«8H  FREG. 

ITITLE(  10)=8h 

ITITLEI 1 1 )=8H     R.H. 

ITITLE(12)=8HMARKS 

CALL  DRAW(  15C,EE,V,0,0,l<H 

END 

END 


,ITITLE,C,C, 4, 0,2,2,9, 8,1, LAST) 


RESCNANT  FREQ 


DAMPING  FACTCR   COEFFICIENT 


VOLUNE  CCNC./ 


1 

3, 

.8665846E+05 

1 *0 1 73467E- 

-01 

4.5119378E- 

■14 

9.50CCCC0E- 

■12 

2 

3 

.7511642E+05 

1  .0074£7i4E- 

-01 

3.  1393589E- 

13 

6.610CCC0E- 

■11 

3 

3, 

.6424348E+C5 

9.9797C43E- 

-02 

5.8275239E- 

■13 

1.227CCCCE- 

■1C 

4 

3, 

.5398310E+05 

9.8883C59E- 

-02 

8.515689CE- 

•13 

1.793CCCCE- 

-1C 

5 

3. 

.UU28493E+05 

9.80C2^92E- 

-C2 

1. 12C3854E- 

•12 

2.3590C00E- 

-1C 

6 

3 

.3510400E+C5 

9.7153255E- 

-02 

1.3S92019E- 

•12 

2.925CCC0E- 

-10 

7 

3. 

.264C000E+05 

9.6333^4CE- 

-C2 

1.658C184E- 

•12 

3.491CCC0E- 

-10 

8 

3, 

.  1813671E+C5 

9.5541  <C3E- 

-02 

1.92683149E- 

•12 

4.057CCC0E- 

-1C 

9 

3, 

.  10281 48E  +  05 

9.47752l^3E- 

-02 

2.  1956514E- 

•12 

4.6230CC0E- 

-1C 

10 

3, 

.0280482E+C5 

9.4033789E- 

-02 

2.46M4680E- 

•12 

5.1890CCCE- 

•1C 

11 

2, 

.9568000E+C5 

9.3315586E- 

-02 

2.7332845E- 

•12 

5.755CCC0E- 

-10 

12 

2. 

.8888276E+05 

9.2619283E- 

-02 

3.CC21G1CE- 

•12 

6.321CCC0E- 

-1C 

13 

2, 

.8239101E+G5 

9.19144C22E- 

-02 

3.27C9175E- 

■12 

6.8870CCOE- 

-1C 

14 

2, 

.7618462E+C5 

9.1288»433E- 

-02 

3.5397340E- 

■12 

7.453CCCCE- 

-1C 

15 

2, 

,702^516E+C5 

9.065U2CE- 

-02 

3.8C855C5E- 

•12 

8.019CCCCE- 

■1C 

16 

2, 

►6455579E+C5 

9.003246CE- 

-02 

4.077367CE- 

■12 

8.5850CC0E- 

-1C 

17 

2 

.591C103E+C5 

8.9430e9CE- 

-02 

4.31461835E- 

•12 

9.1510CC0E- 

■1C 

18 

2. 

.5386667E+C5 

8.88«449C9E- 

-02 

4.615CCCCE- 

•12 

9.717CCCCE- 

•10 

19 

2, 

.4883960E+C5 

8.8274568E- 

-02 

4.8638165E- 

•12 

1.O283CC0E- 

•C9 

20 

2. 

,i4U00777E  +  C5 

8.7718S64E- 

-02 

5.  1526330E- 

•12 

1.0849CC0E- 

-C9 

21 

2, 

.3936000E+C5 

8  .7  1774i43E- 

-02 

5.421HH95E- 

•12 

1.1415CCCE- 

-C9 

22 

2. 

.3488598E+C5 

8.66149289E- 

-02 

5.69C266CE- 

•12 

1.1981CCCE- 

-09 

23 

2, 

.3057615E+C5 

8.6134225E- 

-02 

5.959C826E- 

■12 

1.25U7CCCE- 

-C9 

24 

2 

.2642162E+C5 

8.5631407E- 

-02 

6.2278991E- 

•12 

1.3113CCCE- 

-C9 

25 

2. 

.2241416E+05 

8.51  I401j27E- 

-02 

6.U967156E- 

■12 

1.3679CCCE- 

-C9 

26 

,  1854609E  +  C5 

8.4660EC2E- 

-02 

6.7d55321E- 

■12 

1.4245CCCE- 

-C9 

27 

2. 

.1481026E+C5 

8.4192C79E- 

-02 

7.0343486E- 

■12 

1.4811CC0E- 

-C9 

28 

2, 

.  1  12C000E  +  C5 

8.3733E31E- 

-02 

7.3C31651E- 

■12 

1.5377C0CE- 

-C9 

29 

.077C909E+C5 

8.3285d52E- 

-02 

7.5719816E- 

•12 

1.59U3CC0E- 

-C9 

30 

2, 

.0U33171E+C5 

8.2847162E- 

-02 

7.814C7981E- 

■12 

1.65C9CCCE- 

-C9 

31 

2, 

.01062UOE+C5 

8.2417997E- 

-02 

8.  1C96146E- 

■12 

1.7075CCCE- 

-C9 

32 

1, 

.9789606E+C5 

8.1997E16E- 

-02 

8.  378M31 1 E- 

-12 

1.7641CCGE- 

-C9 

33 

1, 

.9U82791E+05 

8.1586292E- 

-02 

8.6M72U76E- 

•12 

1.8207CCCE- 

-C9 

34 

1, 

.91853U4E+C5 

8.1183117E- 

-02 

8.916C641E- 

■12 

1.8773CCCE- 

-C9 

35 

1 

.8896842E+05 

8.0787999E- 

-02 

9.  1EM88C7E- 

■12 

1.9339CC0E- 

-C9 

36 

1 

.8616889E+05 

8.04C0659E- 

-02 

9.4536972E- 

■12 

1.9905CCCE- 

-C9 

37 

1, 

.8345109E+C5 

8.0020E31E- 

-02 

9.7225137E- 

■12 

2.0471CCCE- 

•C9 

38 

1 

.8G81151E+05 

7.9648264E- 

-02 

9.99133C2E- 

■12 

2.1C37CCCE- 

-C9 

39 

1, 

.782*4681E  +  C5 

7.9282717E- 

-02 

1.C26C147E- 

•11 

2.16C3CC0E- 

•C9 

40 

1, 

.7575385E+C5 

7.8923962E- 

-02 

1.0528963E- 

■11 

2.2169C0CE- 

-C9 

41 

1, 

•7332966E+C5 

7.8571779E- 

-02 

1.0797780E- 

■11 

2.2735CCCE- 

-C9 

42 

1. 

.7097143E+05 

7.822596CE- 

-02 

1.  1C66596E- 

-11 

2.33C1CCCE- 

-C9 

43 

1, 

.6867651E+C5 

7.7886205E- 

-02 

1.  1486C50E- 

-11 

2.4184170E- 

-C9 

44 

1, 

.6644238E+05 

7.7552625E- 

-02 

1. 1C35669E- 

-11 

2.3235681E- 

-C9 

45 

1, 

.6U26667E+05 

7.7224737E- 

-02 

1.060853CE- 

•11 

2.2336529E- 

-C9 

46 

1, 

.621471CE+05 

7.69C2i»67E- 

-02 

1.02C3153E- 

-11 

2.1482998E- 

-C9 

47 

1 

.6CG8153E+05 

7.65856M8E- 

-02 

9.8181697E- 

-12 

2.0672*<C6E- 

-CS 
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48  1.5806792E+05 

49  1.5610435E+C5 

50  1.5U18896E+C5 

51  1.5232000E+C5 

52  1.5C495B1E+05 

53  1.U871U79E+C5 
5U  1  .46975UUE+C5 

55  1  .4527630E  +  05 

56  1  .4361600E  +  05 

57  1.4199322E+C5 

58  1  .14CU0670E  +  C5 

59  1.3885525E+C5 

60  1 .3735770E+05 

61  1.3585297E+C5 

62  1  .344Q000E+05 

63  1  .3297778E  +  C5 

64  1.315E534E+05 

65  1.3022176E+C5 

66  1  .2883615E+05 

67  1.2757766E+05 

68  1.2629548E+C5 

69  1.2503881E+C5 

70  1  .238C690E  +  05 

71  1.2259902E+C5 

72  1.2141449E+C5 

73  1.2025263E+G5 

74  1.1911280E+C5 

75  1 . 1799437E+C5 

76  1.1689674E+05 

77  1 .1581935E+05 

78  1  .  1476164E  +  05 

79  1  .  1372308E  +  C5 

80  1.1270314E+05 

81  1.1  170133E+C5 

82  1.1C71718E+05 

83  1.0975022E+C5 

84  1.03800C0E+C5 
*85  1 .0786609E+C5 

86  1 .069U809E+C5 

87  1.C604557E+05 

88  1.0515816E+C5 

89  1.0U2851J8E  +  C5 

90  1 .0342716E+05 

91  1.C256286E+05 

92  1.0175223E+C5 

93  1  .0C934914E  +  C5 

94  1 .0C13068E+C5 

95  9.9339130L+C4 

96  9.8560000E+C4 

97  9.7792996E+04 

98  9.7C37838E+04 

99  9.6294253E+C4 

100  9.5561977E+C4 

101  9.4840755E+04 

102  9.U13C337E+C4 

103  9.3M30483E+C4 

104  9.27»40959E  +  CU 

105  9.2061538E+04 

106  9. 13920G0E+04 

107  9.0732130E+04 

108  9.0081720E+C4 

109  8.9440569E+04 

110  8.8808481E+04 

111  8.8185263E+04 

112  8.7S70732E+04 

113  8.6964706E+04 

114  8.6367010E+04 

115  8.5777474E+04 


7.6274120E-02 

7.S96773CE-02 

7.S66633CE-02 

7.5369781E-02 

7.5077S^5E-C2 

7.4790£91jE-C2 

7.45C7SC2E-02 

7.42291JU9E-02 

7.3955219E-02 

7.36851C2E-02 

7.3418S9CE-02 

7.3156780E-02 

7.2898372E-C2 

7.2643671E-C2 

7.2392585E-02 

7.2145C24E-02 

7.1900SC2E-02 

7.1660136E-C2 

7.14226U6E-02 

7.1  188255E-02 

7.0957187E-02 

7.0729C71E-C2 

7.05C3c,35E-02 

7.0281713E-02 

7.OO62539E-02 

6.98457149E-02 

6.9631E82E-C2 

6.9i420d78E-02 

6.9212C79E-02 

6.9006C31E-02 

6.88021J77E-02 

6.8601 267E-02 

6.8U02<H*8E-02 

6.8206272E-02 

6.8012189E-02 

6.7820Z55E-02- 

6.7630722E-C2 

6.7443248E-02 

6.7257E89E-02 

6.7074d03E-02 

6.6893352E-02 

6.6714C9UE-02 

6.6536792E-02 

6.6361M10E-02 

6.6187CC9E-02 

6.6016256E-02 

6.5846l»UE-02 

6.5678Z56E-02 

6.5512C43E-02 

6.53U7i4ii6E-02 

6.5184533E-02 

6.5023275E-02 

6.4863£42E-C2 

6.47C56C6E-02 

6.4549129E-02 

6.439U21UE-02 

6.U214  0ECUE-02 

6.4088E84E-02 

6.3938^28E-02 

6.3789412E-02 

6.36»4l611E-02 

6.3495*C2E-02 

6.3350763E-02 

6.320727CE-C2 

6.30651C3E-02 

6.292424CE-02 

6.2784659E-02 

6.2646241E-02 


9.U522135E-12 

9.  1C4U112E-12 

8.7732747E-12 

8. U5e  19U2E-12 

8.  1579319E-12 

7.8717153E-12 

7.59E7326E-12 

7.3382277E-12 

7.089H957E-12 

6.8518791E-12 

6.62U7638E-12 

6.UC7576CE-12 

6. 1997792E-12 

6.0C08712E-12 

5.81C3818E-12 

5.6278703E-12 

5.4529233E-12 

5.2E5153CE-12 

5.  1241951E-12 

4.9697C73E-12 

M.8212678E-12 

4.67E6737E-12 

4.5419298E-12 

U.U1C2977E-12 

4.2836941E-12 

4.  16189C3E-12 

4.0W6610E-12 

3.9217933E-12 

3.823C864E-12 

3.7183502E-12 

3.617MC51E-12 

3.52CC811E-12 

3.4262173E-12 

3.3356612E-12 

3.2U82685E-12 

3.  1639C22E-12 

3.082ti323E-12 

3.0C37358E-12 

2.9276955E-12 

2.8542C03E-12 

2. 7831447E-12 

2.7144282E-12 

2.6U79554E-12 

2.5836355E-12 

2.5213819E-12 

2.461  1  12UE-12 

2.4C27486E-12 

2.3462157E-12 

2.2914424E-12 

2.2382610E-12 

2. 1869C64E-12 

2. 137C17CE-12 

2.0886335E-12 

2.0U16997E-12 

1.9961616E-12 

1.9519678E-12 

1.9C9C690E-12 

1.8674 180E-12 

1.8269700E-12 

1.7876817E-12 

1.7495118E-12 

1.7124209E-12 

1.6763712E-12 

1.6413262E-12 

1.6C7251UE-12 

1.5741132E-12 

1. 54167986-12 

1.51052C56-12 


1.9902C87E-C9 

1.9169569E-C9 

1.8472564E-C9 

1.78C8943E-C9 

1.7176733E-C9 

1 .6574C97E-C9 

1 .5999225E-C9 

1.545C825E-CS 

1.49271146-C9 

1.44266C6E-C9 

1  .39486C9E-C9 

1.3U9131UE-C9 

1.3053793E-C9 

1.263US87E-C9 

1.223Z9C7E-C9 

1.  18U962J4E-C9 

1.1481269E-C9 

1.1128C2UE-C9 

1.C789123E-C9 

1.0U63eU5E-C9 

1.0151513E-C9 

9.851M876E-1C 

9.5631699E-1C 

9.28599141E-1C 

9.019«4271E-1C 

8.7629660E-1C 

8.5161366E-1C 

8.278U9C9E-1C 

8.0U96C57E-10 

7.829C81CE-1C 

7.6165385E-1C 

7.U1 162CUE-10 

7.2139877E-1C 

7.0233196E-1C 

6.8393120E-1C 

6.6616765E-1C 

6.U901397E-1C 

6.324UU21E-1C 

6.  16M337UE-1C 

6.0095915E-1C 

5.85S982CE-1C 

5.7152977E-1G 

5.5753275E-1C 

5.U3991C2E-1C 

5.3088338E-1C 

5.  18193148E-1C 

5.05901483E-10 

U.9U00168E-1C 

U.82U6SC4E-1C 

U.7129260E-1C 

U.60U5871E-1C 

U.4995^36E-1C 

U. 397671 1E-1C 

U.298e5C7E-lC 

U.2029692E-1C 

U.109S179E-10 

U.0195933E-1C 

3.9318962E-1C 

3.8U67318E-1C 

3.76U0C92E-10 

3.6836U17E-1C 

3.6055459E-1C 

3.5296422E-1C 

3.U558542E-1C 

3.38141C86E-1C 

3.31U335UE-1C 

3.2U6U672E-1C 

3.180U394E-10 
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116 
117 
118 
1  19 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
13C 
131 
132 
133 
134 
135 
136 
137 
138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
.149 
150 
151 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 

I  I  I 

172 
173 
174 
175 


7 
7 

7 
7 
7 
7 
7 
7 
7 
7 
7 


7 
7 
7 

7 


5  19^92 
462222 
405618 
349767 
294653 
.2U0262 
.186^79 
133592 
081286 
C29648 
7.978666 
7.928328 
878620 
829532 
.781052 
733169 
685871 
639148 
592990 
547387 
502328 
457804 
413805 
7.370322 
7.327346 
284869 
2^2881 
2C1375 
160341 
7.1 19773 
7.079662 
7.GH0C00 
7.CC0779 
6.961994 
6.923636 
6.885698 
6.848174 
6.811056 
6.774339 
6.738016 
6.7C2080 
6.666525 
6.631345 
6.596535 
6.562C88 
6.528000 
6.494263 
6.460874 
6.427326 
6.395114 
6.362734 
6.330680 
6.298947 
6.267531 
6.236426 
o  •  ^Uaoz'y 
6.175135 
6.144938 
6.1 15C36 
6.085423 


2E  +  C4 
26  +  04 
7E+C4 
4E  +  C4 
5E  +  C4 
3E  +  C4 

ec+C4 

2E  +  C4 
2E  +  C4 
6E  +  C4 
7E  +  C4 
1E  +  C4 
7E  +  C4 
7E  +  04 
6E  +  C4 
2E  +  C4 
6E  +  04 
9E  +  C4 
9E  +  C4 
4E  +  C4 
4E  +  C4 
2E  +  C4 
3E  +  C4 
6E  +  C4 
9E  +  04 
6E  +  04 
8E  +  C4 
4E  +  C4 
9E  +  04 
4E  +  C4 
GE  +  C4 
CE  +  C4 
9E  +  04 
5E  +  04 
4E  +  C4 
6E  +  04 
4E  +  04 
9E  +  04 
6E  +  C4 
1E  +  04 
CE+04 
2E  +  04 
6E  +  C4 
4E  +  C4 
8E  +  C4 
0E+C4 
6E  +  04 
0E  +  04 
1E  +  04 
5E  +  04 
2E  +  G4 
1E  +  04 
4E  +  04 
2E  +  04 

8E  +  C4 

/  —  .  /« «. 

Oc  +  OM 

1E  +  04 
9E  +  C4 
5E  +  C4 
7E  +  04 


6.25C9 
6.2373 
6.2238 
2105 
1973 
1841 
1711 
1582 
1455 
1328 
1202 
1078 
.0954 
.0832 
.0710 
6.0589 
6.0470 
6.0351 
6.0234 
6.0117 
6.00C1 
5.9886 
5.9772 
5.9659 
5.9547 
5.9435 
5.9325 
5.9215 
5.91C6 
.8998 
.8891 
8785 
8679 
.8574 
5.8U70 
5.8367 
5.8264 
5.8162 
5.8061 
5.7961 
5.7861 
5.7762 
5.7664 
5.7566 
5.7469 
5.7373 
5.7277 
5.7182 
5.7083 
5.6994 
5.6901 
5.68C8 
5.6717 
"5.6625 
5.6535 
-5.6445 
5.6355 
5.6266 
5.6178 
5.6090 


265E-02 
412E-02 
76ME-C2 
3C1E-02 
CC5E-02 
E59E-02 
6UUE-02 
945E-02 
143E-02 
l^24E-02 
771E-02 
168E-C2 
599E-02 
C51E-02 
5C8E-C2 
956E-02 
281E-C2 
769E-02 
1C6E-02 
379E-02 
575E-C2 
682E-02 
6E7E-C2 
577E-C2 
341E-02 
967E-02 
^43E-C2 
758E-02 
S01E-C2 
861E-C2 
627E-02 
189E-02 
536E-02 
658E-C2 
5i*6E-02 
189E-02 
578E-02 
704E-C2 
558E-C2 
129E-02 
41CE-02 
392E-02 
C65E-0? 
IJ22E-02 
455E-C2 
154E-C2 
513E-02 
523E-02 
176E-C2 
465E-C2 
382E-02 
921E-02 
C73E-02 
E32E-C2 
19CE-C2 
i4  ic-u2 
678E-02 
795E-02 
484E-02 
739E-02 


1.4EC 
1.450 
1.421 
1.393 
1.365 
1.339 
1.312 
1.287 
1.262 
229 
215 
192 
17C 
148 
127 
1C6 
086 
066 
0U7 
028 
C1C 
9.927 
9.752 
9.582 
9.415 
9.252 
9.093 
8.938 
8.786 
8.637 
8.492 
8.35C 
6.211 
8.076 
7.943 
7.813 
7.686 
7.562 
7. 440 
7.321 
7.2C5 
7.09C 
6.979 
6.869 
6.762 
6.658 
6.555 
6.454 
6.356 
6.259 
6.  165 
6.072 
5.981 
5.892 
5.  805 
■5.715 
5.625 
5.553 
5.472 
5.393 


CC59Erl2 
<G77E-12 
2987E-12 
2531E-12 
8456E-12 
1523E-12 
1 50 1 E-  12 
8167E-12 
12C8E-12 
C718E-12 
6200E-12 
7563E-12 
4622E-12 
7206E-12 
5140E-12 
8261E-12 
6414E-12 
9444E-12 
7207E-12 
9560E-12 
6367E-12 
4976E-13 
8242E-13 
2248E-13 
58 126-13 
7794E-13 
7C92E-13 
2647E-13 
2428E-13 
81444E-13 
6736E-13 
7377E-13 
9471E-13 
2153E-13 
4582E-13 
5950E-13 
5472E-13 
2389E-13 
5967E-13 
5494E-13 
C283E-13 
9666E-13 
2999E-13 
9656E-13 
9C32E-13 
0538E-13 
3607E-13 
7687E-13 
2244E-13 
6758E-13 
0727E-13 
3664E-13 
5C96E-13 
4563E-13 
1621E-13 
5637E-13 
6792E-13 
4C78E-13 
7300E-13 
6075E-13 


3.1161 
3.0536 
2.9927 
2.9335 
2.8758 
2 • 8 1 S6 
2.764e 
2.7115 
2.6595 


6088 
5595 
.5113 
4641* 
4186 
2.374C 
2.3304 
2.2879 
2.2464 
2.206C 
2.1664 
2.1279 
2.0902 
2.0534 
2.0175 
1.9824 
1.9481 
9147 
8819 
8499 
8187 
7881 
1.7582 
1.7290 
1.7004 
1.6725 
1.6451 
1.6184 
1.5922 
1.5666 
1.5415 
1.517C 
4930 
4695 
4464 
4239 
4018 
1.38C2 
1.3590 
1.3383 
1.3179 
1.29eC 
1.2785 
1.2594 
1.24C6 
1.2222 
1.2C42 
1.1866 
1.1692 
1.1522 
1.1356 


9C1E-10 
597E-1C 
912E-1C 
298E-1C 
227E-1C 

7C9E-1C 
3C8E-1C 
54CE-1C 
972E-1C 
188E-1C 
787E-1C 
3E2E-1C 
6C4E-1C 
C93E-1C 
5C6E-1C 
511E-1C 
786E-1C 
C25E-1C 
930E-1C 
213E-1C 
599E-1C 
82CE-1C 
618E-1G 
7H6E-1C 
962E-1C 
C36E-10 
744E-10 
868E-1C 
2C1E-1C 
54CE-1C 
691E-1C 
464E-1C 
677E-1C 
153E-1C 
723E-1C 
221E-1C 
486E-1C 
366E-1C 
7C9E-1C 
27CE-1C 
21 1E-1C 
C94E-10 
887E-1C 
465E-1C 
7C2E-10 
479E-1C 
680E-1C 
192E-1G 
9C7E-1C 
718E-10 
522E-1C 
22 1E-1C 
717E-1C 
917E-1C 
72QE-1C 
C66E-1C 
842E-10 
972E-10 
378E-1C 
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N/ELCCITY    SFIFT    CUE    TC    EUBBLES 
IN    MCNEIL    DISTRIBUTION 


FREQUENCY 

RHO  X  COMPRESSIBILITY 

VELOCITY 

1.00 

KC 

1.58C5603E-1 

2 

-.0C26672 

2.00 

KC 

6.4159358E-1 

2 

-.01C8234 

3.0G 

KC 

1  .47969C4E-1 

1 

-.0249692 

4.00 

KC 

2.72C9566E-1 

1 

-.046CB26 

5.00 

KC 

4.5CC6175E-1 

1 

-.0759422 

6.00 

KC 

6.99C7125E-' 

1 

-.  1  1?9544 

7.00 

KC 

1  .0658786E-1 

C 

-.  1798347 

8.00 

KC 

1 .7CE8791E-1 

C 

-.2874469 

9.00 

KC 

2.2668272E-' 

c 

-.3823808 

10.00 

KC 

2.0659377E-1 

c 

-.3M85055 

1  1.00 

KC 

2.0567619E-1 

C 

-.3469582 

12. OC 

KC 

2.031335OE-1 

C 

-.351  1C20 

13.00 

KC 

2.1C27664E-' 

c 

-.3547160 

14.00 

KC 

2.1C49692E-1 

C 

-.355C874 

15.00 

KC 

2.0777052E-1 

c 

-. 35CM899 

16.00 

KC 

2.0126459E-" 

C 

-.3395187 

17.00 

KC 

1  .90188C2E-' 

IC 

-.32C8393 

13.00 

KC 

1  .7371231 E- 

C 

-.2930536 

19.00 

KC 

1  .5C911 19E-' 

IC 

-.2545978 

20.00 

KC 

1  .2C69606E- 

IC 

-.2C36332 

21.00 

KC 

8.1730041E- 

1 1 

-. 1379CC4 

22.00 

KC 

3.2297022E-' 

1 1 

-.0544983 

23.00 

KC 

-2.9915002E- 

1 1 

.05C4841 

24.00 

KC 

-1  .0S192C6E-' 

IC 

. 1826074 

25.00 

KC 

-2.07M24U0E- 

IC 

.3501512 

26.00 

KC 

-3.3431118E-" 

IC 

.561)4686 

27.00 

KC 

-4.86241  15E- 

IC 

.8212C58 

28.00 

KC 

-6.2C75042E-' 

IC 

1.0486149 

29.00 

KC 

-7.18528520- 

IC 

1.2139890 

30.00 

KC 

-7.8912450E- 

IC 

1.3334235 

31.00 

KC 

-8.4C72293E- 

!C 

1.42C7359 

32.00 

KC 

-8.7863755E- 

IC 

1.4849C29 

33.00 

KC 

-9.0640521E-' 

IC 

1.5219023 

'  34.00 

KC 

-9.2647964E-' 

IC 

1.5658830 

35.00 

KC 

-9.4C62358E-' 

IC 

1.5898263 

36.00 

KC 

-9.5C13636E- 

IC 

1.6059204 

37.00 

KC 

-9.56CC018E-' 

IC 

1.6158575 

38.00 

KC 

-9.5896551E-" 

IC 

1.6208778 

39.00 

KC 

-9.5961931E-' 

IC 

1.6219846 

40.00 

KC 

-9.5841798E- 

iC 

1.61995C8 

41.00 

KC 

-9.5572699E-^ 

IC 

1.6153950 

42.00 

KC 

-9.5183848E-' 

IC 

1.6088120 

43.00 

KC 

-9.4698572E- 

IC 

1.6CC5967 

44.00 

KC 

-9.4135676E- 

IC 

1.591C674 

45.00 

KC 

-9.351C402E-' 

IC 

1.58C4825 

46.00 

KC 

-9.2835111E-- 

IC 

1.569C510 

47.00 

KC 

-9.21  19815E- 

IC 

1.5569426 

48.00 

KC 

-9.1372639E- 

IC 

1.5442949 

49.00 

KC 

-9.06CC184E-' 

IC 

1.5312196 

50.00 

KC 

-8.9807696E-" 

IC 

1.5178C55 

51  .00 

KC 

-8.8999037E- 

IC 

1.5C41  181 

52.00 

KC 

-8.8176999E- 

IC 

1.49C2046 

53.00 

KC 

-8.7344152E-" 

IC 

1.4761086 

54.00 

KC 

-8.65C2427E-' 

IC 

1.4618627 

55.00 

KC 

-8.5651873E- 

IC 

1.4474678 

56.00 

KC 

-8.4792771E-' 

IC 

1.4329287 

57.00 

KC 

-8.3925645E- 

IC 

1.4182542 

58.00 

KC 

-8.30M8841E- 

IC 

1.4C3M163 

59.00 

KC 

-8.21631 18E- 

IC 

1.3884280 

60.00 

KC 

-8.1269672E- 

IC 

1.3723C94 

61.00 

KC 

-3.0375212E- 

IC 

1.3581741 

62.00 

KC 

-7.9495245E- 

IC 

1.3432845 

63.00 

KC 

-7.8652497E- 

IC 

1.229C251 

64.00 

KC 

-7.7861964E- 

IC 

1.3156495 

65.00 

KC 

-7.7129255E- 

IC 

1.3032527 

66.00 

KC 

-7.6453344E- 

IC 

1.2918171 

67.00 

KC 

-7.5829973E- 

IC 

1.2812706 

M/S 

v/s 

N/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
f*/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
t*/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
M/S 
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FREQUENCY 

RHO  X  CCMPRESSII 

68. OC 

KC 

-7.5253968E-1C 

69.00 

KC 

-7.472C296E-1C 

70.00 

KC 

-7.4224420E-1C 

71.00 

KC 

-7.3762367E-1C 

72.00 

KC 

-7.333C692E-1C 

73.00 

KC 

-7.29264  1 1 6-1 C 

74.00 

KC 

-7.2546934E-1C 

75.00 

KC 

-7.219C0C0E-1C 

76.00 

KC 

-7.1853630E-1C 

77.00 

KC 

-7.1536079E-1C 

78.00 

KC 

-7.12358C5E-1C 

79.00 

KC 

-7.0951439E-1C 

80.00 

KC 

-7.0681755E-1C 

81.00 

KC 

-7.0425657E-1C 

82.00 

KC 

-7.0182160E-1C 

83.00 

KC 

-6.9950375E-1C 

84.00 

KC 

-6.9729491E-K 

85.00 

KC 

-6.9518782E-1C 

86.00 

KC 

-6.9217579E-1C 

87.00 

KC 

-6.9125275E-1C 

88.00 

KC 

-6.8941313E-1C 

89.00 

KC 

-6.8765182E-1C 

90.00 

KC 

-6.8596414E-1C 

91.00 

KC 

-6.8434574E-1C 

92.00 

KC 

-6.8279266E-1C 

93.00 

KC 

-6.812C120E-1C 

94.00 

KC 

-6.7986794E-1C 

95.00 

KC 

-6.7848972E-1C 

96.00 

KC 

-6.7716359E-1C 

97.00 

KC 

-6.7588683E-1C 

98.00 

KC 

-6.7465687E-1C 

99.00 

KC 

-6.7347134E-1C 

100.00 

KC 

-6.7232802E-1C 

.101.00 

KC 

-6.7122482E-1C 

102.00 

KC 

-6.7C15982E-1C 

103.00 

KC 

-6.69131  17E-1C 

104.00 

KC 

-6.6813717E-1C 

105.00 

KC 

-6.6717622E-1C 

106.00 

KC 

-6.6624680E-1C 

107.00 

KC 

9  -6.6524749E-1C 

108.00 

KC 

-6.6447695E-1C 

109.00 

KC 

-6.6363392E-1C 

110.00 

KC 

-6.6281720E-1C 

111.00 

KC 

-6.62C2568E-1C 

112.00 

KC 

-6.6125829E-1C 

113.00 

KC 

-6.6051402E-1C 

114.00 

KC 

-6.5979192E-1C 

115.00 

KC 

-6.59C91C9E-1C 

116.00 

KC 

-6.5841067E-1C 

117.00 

KC 

-6.5774985E-1C 

118.00 

KC 

-6.571C787E-1C 

119.00 

KC 

-6.5648399E-1C 

120.00 

KC 

-6.5587751E-1C 

121.00 

KC 

-6.5528779E-1C 

122.00 

KC 

-6.5471417E-1C 

123.00 

KC 

-6.5415608E-1C 

124.00 

KC 

-6.5361293E-1C 

125.00 

KC 

-6.5308418E-1C 

126.00 

KC 

-6.5256932E-1C 

127.00 

KC 

-6.52C6785E-1C 

128.00 

KC 

-6.5157929E-1C 

129.00 

KC 

-6.5110321E-1C 

130.00 

KC 

-6.5063915E-1C 

131.00 

KC 

-6.5018672E-10 

132.00 

KC 

-6.4974552E-1C 

133.00 

KC 

-6.4931517E-1C 

134.00 

KC 

-6.4889531E-1C 

VELOCITY    SHIFT 


1.2715257 

M/S 

1.2624971 

M/S 

1.2541C61 

M/S 

1.2462915 

M/S 

1.2389888 

m/s 

1.2221497 

M/S 

1.2257303 

M/S 

1.2196922 

M/S 

1.2140022 

M/S 

1.2086205 

M/S 

1.2C35512 

M/S 

1.  1987410 

M/S 

1. 1941791 

M/S 

1. 1898472 

M/S 

1. 1857284 

M/S 

1. 1818C77 

M/S 

1.  178C716 

M/S 

1. 1745074 

M/S 

1. 1711C42 

M/S 

1. 1678515 

M/S 

1. 1647399 

M/S 

1. 1617607 

M/S 

1. 1589C61 

M/S 

1.  1561688 

M/S 

1.  1535419 

M/S 

1. 151C192 

M/S 

1.  1485951 

M/S 

1.  1462640 

M/S 

1.  1440210 

M/S 

1. 1418615 

M/S 

1.  1297812 

M/S 

1. 1277761 

M/S 

1.  1358423 

M/S 

1. 1329765 

M/S 

1.  1221752 

M/S 

1. 1204254 

M/S 

1. 1287543 

M/S 

1.1271290 

M/S 

1. 1255571 

M/S 

1. 124C261 

M/S 

1. 1225627 

M/S 

1.121  1379 

M/S 

1. 1197566 

M/S 

1.  1  184179 

M/S 

1. 11712C1 

M/S 

1.  1158613 

M/S 

1.  1  146400 

M/S 

1. 1124548 

M/S 

1. 1122040 

M/S 

1.111  1864 

M/S 

1. 1  101007 

M/S 

1. 1090455 

M/S 

1. 108C198 

M/S 

1. 1C70224 

M/S 

1. 1060523 

M/S 

1. 1051084 

M/S 

1. 1C41898 

M/S 

1.  1032956 

M/S 

1. 1024249 

M/S 

1. 1015767 

M/S 

1. 1CC75C5 

M/S 

1.0999453 

M/S 

1.0991605 

M/S 

1.0983954 

M/S 

1.0976492 

M/S 

1.0969214 

M/S 

1.0962113 

M/S 
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FREQUENCY 


135.00 
136.00 
137.00 
138.00 
139.00 
140.00 
141.00 
1U2.CC 
143.00 
144.00 
145.00 
146. OC 
147.00 
148.00 
149.00 
150.00 
151.00 
152.00 
153.00 
154.00 
155.00 
156.00 
157.00 
158.00 
159.00 
160.00 
161.00 
162.00 
163.00 
164.00 
165.00 
166.00 
167.00 
168.00 
'169.00 
170.00 
171.00 
172.00 
173.00 
174.00 
175.00 
176.00 
177.00 
178.00 
179.00 
180.00 
181.00 
182.00 
183.00 
184.00 
135.00 
186.00 
187.00 
188.00 
189.00 
190. OC 
191.00 
192.00 
193.00 
194.00 
195.00 
196.00 
197.00 
198.00 
199.00 
200.00 


KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 

KC 
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APPENDIX  VIII 
ENVELOPE  DETECTOR  AND  ECHO  COUNTER 

This  system  starts  a  Beckman  7360  one  megacycle  counter  at  the 
start  of  the  main  bang  and  stops  it  at  the  beginning  of  a  preselected 
echo.  Figure  (IX-2)  is  a  circuit  diagram  of  the  system. 

The  system  operation  is  as  follows.  The  envelope  detector 
passes  the  envelope  of  the  echo  train  but  does  not  pass  the  high 
frequency  signal.  The  "and"  circuit  passes  the  echo  train  only 
between  the  main  bang  and  the  preselected  echo  for  which  we  are 
measuring  travel  time.  The  amplifier  and  clipper  convert  the  echo 
train  into  a  rectangular  pulse  train  with  steep  leading  and  trailing 
edges.  The  echo  counter  counts  echos  and  sends  a  stop  pulse  to  the 
one  megacycle  counter  when  the  echo  counter  reaches  the  preselected 
count.  At  this  time  flip  flop  one  is  reset  by  this  stop  signal. 
The  flip  flop  output  turns  off  the  clipper-amplifier.  The  next 
bang  starts  the  one  megacycle  counter  and  the  flip  flo  which  in 
turn  turns  the  amplifier-clipper  back  on.  The  process  then 
repeats  itself. 

The  system  as  shown  is  plagued  by  two  sources  of  error.  First, 
as  frequency  increases  the  rise  time  of  the  detected  wave  pulse 
decreases  thus  affecting  the  counter  reading.  Second^  as  attenuation 
increases  the  rise  time  of  the  detected  wave  pulse  increases 
affecting  the  counter  reading.  The  combination  of  these  effects 
caused  velocity  variation  of  five  percent  from  15  KC  to  100  KG  over 
two  round  trips  in  the  laboratory  "non-dispersive"  water. 
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By  using  a  more  elaborate  clipping  network  with  high  sensitivity 
perhaps  the  above  effects  could  be  minimized.  The  wave  shape  of 
the  detected  pulse  train  limited  the  gain  that  could  be  used  in  the 
clipper  as  shown  in  figure  (IX-l).  This  problem  might  be  solved  by 
using  several  clippers  with  different  sensitivity.  The  most  sensi- 
tive clipper  would  be  turned  on  only  when  the  last  echo  to  be  counted 
is  received.  This  would  enable  later  echos  to  be  received. 
Receiving  the  16th  echo  may  decrease  the  variation  of  0.6  percent 
and  increasing  the  gain  of  this  last  echo  may  decrease  the  vari- 
ation further.  Figure  (IX-3)  is  a  possible  block  diagram  of  the 
proposed  improved  system.  Time  did  not  permit  construction  of 
the  improved  system. 
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APPENDIX  IX 
TIME  MARKER 

This  system  presents  a  display  on  a  dual  trace  oscilloscope 
similar  to  a  radar  "A"  scope.  A  selected  echo  is  superimposed  on 
a  step  voltage.  The  step  can  be  moved  to  any  point  on  the  echo. 
This  step  stops  a  one  megacycle  counter  which  was  started  by  the 
main  bang. 

The  heart  of  the  system  is  a  DuMont  404  B  pulse  generator 
with  delay.  Delays  from  two  microseconds  to  12.5  milliseconds  can 
be  obtained  in  four  ranges.  With  the  dual  trace  oscilloscope 
synchronized  to  the  start  of  the  delayed  pulse  any  echo  up  to  the 
11th  echo  can  be  displayed  in  detail.  The  pulse  width  can  be  finely 
controlled  with  a  four  inch  diameter  knob  providing  100  microseconds 
change  for  a  3&0  degree  turn.  The  rise  and  fall  times  of  the  pulse 
are  14  nanoseconds  maximum.  Pulse  widths  from  .05  to  105  micro- 
seconds in  three  ranges  are  available. 

The  Beckman  7360  one  megacycle  counter  is  started  by  the  gating 
voltage  output  of  the  tone  burst  generator.  This  gating  voltage 
output  has  an  amplitude  of  24  volts,  and  width  determined  by  the 
main  bang  width.  The  counter  is  stopped  by  the  trailing  edge  of  the 
pulse  generator  ouput.  The  counter  can  be  started  and  stopped  by 
a  minimum  plus  or  minus  one  volt  pulse. 

Overall  error  due  to  instabilities  in  the  counter,  pulse 
generator,  and  oscilloscope  are  indicated  in  figure  (IX-l).  A 
mean  deviation  of  about  .011  percent  is  obtained  from  Figure 
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(iX-l)  when  using  the  seventh  echo  for  measurements  and  .0094 
percent  when  using  the  eleventh  echo.  When  taking  actual  measure- 
ments the  error  can  be  reduced  by  visual  integration  of  the  counter 
readout. 
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APPENDIX  X 

Extracts  from  N.E.L.  Tower  log  together  with  temperature  and 

oxygen  data. 

10  February  1965 

Tides:  High  1746  at  3.2  ffcrj  Lows  1119  at  0.2  ft.,  2202  at  2.4  ft, 

Wave  Height:  3s  ft.;  Water  Clarity  15  ft.,  Slicks,  None. 

Wind:  Speed  0-6  kts.;  Direction  from  100°  to  140°. 

Depth  10  ft.,  Time  1330  -  1614. 

Temp.:  13.2  +  0.1°C;  Oxygen  Content  96  +  3%   saturation. 

Depth  10-45  ft.;  Time  1835  -  1855. 

Temp.  13.25  ±   0.05°C;  Oxygen  Content  101  +  %   saturation 

Depth  45-10  ft.;  Time  1918  -  1958 

Temp.  13.3  ±  0.1°C;  Oxygen  Content  98  ±  l&   saturation 

Depth  10  ft.;  Time  2141  -  2353. 

Temp.  13.4  +  0.1°C;  Oxygen  Content  100.5  ±  1.5%   saturation 

12  February  1965 

Water  Clarity,  20  ft.;  Wind,  Speed  0-3  knots;  Direction  from  095° 

Depth  10-45  ft.;  Time  1156  -  1218 

Temp.  13.2  +  0.1OC;  Oxygen  Content  99.5  ±  3.5$  saturation 

Depth  45-10  ft.;  Time  1224  -  1247- 

Temp.  13-15  ±  0.05°C.;  Oxygen  Content  95  ±  2%   saturation 
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15  February  1965 

Tides;  High  0805  at  7.2  ft.;  Low  1512  at-1.9  ft.;  High  2119  at  5.0  ft, 

Wave  Height  4  ft.;  Water  Clarity  20  ft.;  Slicks  -  None 

Wind;  Speed  2  kts.;  Direction  from  020° 

Depth  10-45  ft.;  Time  0925  -  0942 

Temp.  12.85  +  0.05°C.;  Oxygen  Content  98  +  2%   saturation 

Depth  45-10  ft., |  Time  0950  -  1008 

Temp.  Isothermal  at  12.9°C.;  Oxygen  Content  94.5  ±  3.5$  saturation 

Depth  25  ft.;  Time  1608  -  1639 

Temp.  1305  ±   0.5°C.;  Oxygen  Content  94+1$  saturation 

Oxygen  data  is  not  available  after  this  time  since  the  oxygen  probe 
lead  split  and  water  entered  the  system. 

16  February  1965 

Tides;  Low  0252  at  0.8  ft.;  High  0849  at  7.0  ft.;  Low  1549  at-1.7  ft. 
Wave  Height  2g  ft.;  Water  Clarity  12  ft.;  Slicks ,  Some. 
Wind;  Speed  5  knots;  Direction  from  095° • 
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APPENDIX  XI 
Effect  of  Scatterer  Position  on  Signal 
Experimentation  was  carried  out  to  determine  the  effect  on 
backscattered  sound  from  a  1.5  inch  styrofoam  sphere  when  it  was 
placed  in  a  number  of  positions  within  the  region  between  the 
transducer  and  reflector. 

The  sphere  was  suspended  in  the  space  at  the  positions  shown 
in  Figure  XI-4  and  the  set  of  accompaning  photos  were  taken  to 
show  the  backscatter.  A  partial  analysis  follows: 
Pulse  length  =75.0  cm  -  2.46  ft  =0.5  msec 
Cycles  per  pulse  (at  49.05  KG)  =  24.5 
The  region  is  broken  into  7  approximately  equal  intervals  as 
shown  in  Figure  XI-4.  The  distance  from  the  transducer  to  the 
reflector  is  76.4  cm  or  i=*s  pulse  length. 
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Figure  XI-1 


The  first  backscatter  from  sphere  will  appear  at  t  =  4/7(25)  =  14.3 
cycles  having  travelled  to  position  1  and  back  to  the  transducer 
but  it  will  be  blanked  by  the  transmitted  pulse  for  3/7(25)  =  10.7 
cycles.  Although  the  reflected  pulse  exists  for  a  full  25  cycles 
it  will  be  seen  only  after  the  source  pulse  has  passed  by  the 
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sphere,  that  is  it  will  be  seen  for  only  t  =  4/7(25)  =  14»3 
cycles.  This  is  not  completely  seen  in  photo  1  because  of  the 
extra  wide  blanking  pulse  applied  to  the  main  bang  interval. 
There  is  then  a  break  of  3/7(25)  =  10.7  cycles  =  0.214  msec  before 
the  reflection  from  the  relfector  can  be  seenj  this  checks  with 
photo  1.  The  reflector  echo  also  lasts  for  25  cycles.  To  be  more 
exact,  however,  any  echo  will  persist  for  a  time  equivalent  to 
the  length  of  the  source  pulse  plus  the  length  of  or  thickness  of 
the  insonified  scattererj  that  is  the  sphere  echo  will  persist  for 
75  +  2.54(15)  =  78.1  cm  =  26  cycles.  A  breakdown  of  photo  1  is 
shown  in  Figure  Xl-2. 
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Analysis  of  Photo  1 
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Position  2; 

At  position  3  (3/7  along  the  axis),  the  situation  is  shown 
in  Figure  XI-3. 
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Figure  XI-3 
Sphere  at  Position  2 

The  scattering  from  the  sphere  will  begin  at  6/7(25)  = 
21.4  cycles  =  0.429  msec.  It  will  persist  for  26  cycles  as 
before  and  will  be  followed  by  a  "no  sound"  interval  of  l/7(25) 
=3.57  cycles  =  0.0715  msec;  this  checks  with  photo  2. 

Position  3' 

In  this  position  (4/7  along  the  axis),  the  backscattered 
sound  begins  at  8/7(25)  =  28.6  cycles  =  0.561  msec.  There  is  no 
blanking  by  the  source  pulse  since  it  terminated  0.0714  msec 
earlier.  The  scattered  sound  continues  until  15/7(0.5)  =1.07 
msec  but  will  be  overlapped  from  15/7(0.5)  =  1.07  msec  or  for  the 
last  0.07  msec  of  the  scattered  pulse. 

The  analysis  can  be  continued  in  a  similar  manner  for  positions 
4  and  5. 
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Photos  6  through  8  show  the  backscatter  for  corresponding 
positions  of  the  sphere. 

Taking  the  average  amplitude  of  the  scatter  from  position 
3  as  a  reference  it  was  found  that  for  the  data  taken  the  total 
deviation  from  this  value  throughout  the  region  was  approximately 
10*  i 
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